Draft version November 30, 2010 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



THE INTERNAL-COLLISION-INDUCED MAGNETIC RECONNECTION AND TURBULENCE (ICMART) 

MODEL OF GAMMA-RAY BURSTS 

Bing Zhang 1 , Huirong Yan 2,3 

Draft version November 30, 2010 

ABSTRACT 

The recent Fermi observation of GRB 080916C shows that the bright photosphere emission asso- 
ciated with a putative fireball is missing, which suggests that the central engine likely launches a 
Poynting-flux-dominatcd outflow. We propose a model of gamma-ray burst (GRB) prompt emission 
in the Poynting-flux-dominated regime, namely, the Internal-Collision-induced MAgnetic Reconnec- 
tion and Turbulence (ICMART) model. It is envisaged that the GRB central engine launches an 
intermittent, magnetically-dominated wind, and that in the GRB emission region, the ejecta is still 
moderately magnetized (e.g. 1 < a < 100). Similar to the internal shock (IS) model, the mini-shells 
interact internally at the radius R\s ~ T 2 cAt. Most of these early collisions, however, have little 
energy dissipation, but serve to distort the ordered magnetic field lines entrained in the ejecta. At a 
certain point, the distortion of magnetic field configuration reaches the critical condition to allow fast 
reconnection seeds to occur, which induce relativistic MHD turbulence in the interaction regions. The 
turbulence further distorts field lines easing additional magnetic reconnections, resulting in a runway 
release of the stored magnetic field energy (an ICMART event). Particles are accelerated either di- 
rectly in the reconnection zone, or stochastically in the turbulent regions, which radiate synchrotron 
photons that power the observed gamma-rays. Each ICMART event corresponds to a broad pulse in 
the GRB light curve, and a GRB is composed of multiple ICMART events. This model retains the 
merits of the IS and other models, but may overcome several difficulties/issues faced by the IS model 
(e.g. low efficiency, fast cooling, electron number excess, Amati/Yonetoku relation inconsistency, and 
missing bright photosphere) . Within this model, the observed GRB variability time scales could have 
two components, one slow component associated with the central engine time history, and another 
fast component associated with relativistic magnetic turbulence in the emission region. The model 
predicts a decrease of gamma-ray polarization degree and E p in each ICMART event (broad pulse) 
during the prompt GRB phase, as well as a moderately magnetized external reverse shock. The model 
may be applied to the GRBs that have time-resolved, featureless Band-function spectra, such as GRB 
080916C and most GRBs detected by Fermi LAT. 

Subject headings: gamma-rays burst: general - gamma-ray bursts: individual: GRB 080916C - mag- 
netic reconnection - turbulence 



1. INTRODUCTION 

The composition of the Gamma-Ray Burst (GRB) 
ejecta has remained a mystery until recently. The un- 
certainty lies in the lack of knowledge about the ratio 
between Poynting flux and matter (baryonic) flux 4 , i.e. 5 



Fb AirTpc 2 Airp'c 2 ' 

where B and p are the magnetic field strength and 
matter density in the rest frame of the central engine 
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4 Strictly speaking, the matter flux is the sum of the bary- 
onic flux and the leptonic flux. It is usually dominated by the 
baryonic flux unless the pair number density is so high that 
N±/Nf > > m v /rn e is satisfied. 

5 The last equation applies to the case that the magnetic field 
lines are perpendicular to the direction of the motion, which is 
the configuration invoked in the ICMART model proposed in this 
paper. If magnetic fields are generated in the internal shocks due 
to plasma instabilities, the comoving field lines would have random 
orientations. In this case, one has a ~ B' 2 /8irp' c 2 . 



(or the lab frame), and B' and p' are the corresponding 
quantities in the rest frame comoving with the ejecta. 
The stand a rd pi ct ure is the "fireball" shock model 
dPaczvnskil 119861: iGoodmanl 119861: IShemi and Piranl 
1990tlRees and Meszarodll99alMeszaros and Reedll993t 
Rees and Meszaroslll994ft . Within such a picture, an ini- 
tially hot fireball composed of photons, electron-positron 
pairs, and a small amount of baryons first converts most 
of its thermal energy into kinetic energy, and then dis- 
sipates the kinetic energy in the internal (or sometimes 
external) shocks to power the observed GRB emission. 
Within such a scenario, the magnetic field is assumed 
not to play a dynamically important role in the ejecta, 
i.e. a <C 1. Such a field can be generated in- situ via 
plasma instabilities (IWeibell 119591 : iMedvedev and Loebl 
119991: iNishikawa et all 120051 120091: (Spitkovskvi 12008ft in 
relativistic shocks, in whic h Fermi-accelerated e l ectrons 
cool v i a synchrotron /jit ter (|Meszaros et al.lll994rTavanil 
119961 : Medvedev 2000) or sy n chrotron self-Compton 
(SSC ) (jMeszaros et all 119941 : Kumar and McMahonl 
2008) radiation. An alternative view is that the GRB 
ejecta carries a dynamically important magnetic field 
component, i.e. a 3> 1. The GRB radiation is powered 
by dissipation of the magnet i c field energy in the ejecta 
(|Usovlll992t lThompsorJfl994 iMeszaros and Rees|[l997b1: 
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Lvutikov and Blandfordl 20031: IVlahakis and Konigll 
20031 lLazarian et al. l i27)03l : lL7uti"kovll2006al) . 



Until recently, it has been difficult to diagnose the com- 
position of GRB ejecta. Regardless of the a values, the 
late time afterglow behavior appears the same, which 
is the emission of the "forwardly" shocked circumburst 
medium. Successful modeling of some afterglow data 
therefore does not shed light onto the composition of the 
ejecta. An important phase is when the ejecta energy 
is transferred to the circumburst medium. This is usu- 
ally accom panied by passing of a reverse shock across 
the e jecta (jSari and Piranl [1995: Z hang and Kobavashil 
2005). The brightness of the reverse shock emission 
is found to be dependen t on the magnetization pa- 
rameter (iFan et al.l l2004at IZhang and Kobavashil 120051: 
iMizuno et al.l l2009at iMimica et all I2009D . However , 
both low-cr and high-cr flo ws can lead to a relatively 
dim reverse shock emission ([Zhang and K obavashl l2005l : 
Din and Fanll2007t ). so that it is difficult to robustly con- 
strain a from the observations of the reverse shock emis- 
sion. Finally, the prompt emission light curves and the 
narrow-band spectra (in the energy bands of the previous 
gamma-ray detectors such as BATSE, Swift, etc) may be 
accommodated within the frameworks of both the low-cr 
(internal shock) and high-cr (magnetic dissipation) mod- 
els, so that they do not carry adequate information to 
diagnose GRB composition. The progress in constrain- 
ing GRB composition therefore has been slow. 

In any case, it has been long speculated that the 
GRB ejecta are somewhat magnetized. The putative 
GRB cen tral engines , eith er a black hole - torus sys- 
tem (e.g. iProga et~all [20 03) or a rapidly spinning neu- 
tron star (e.g. lUsovl 11991 . are very likely magnetized. 
The GRB ejecta is likely entrained with a globally struc- 
tured magnetic field. In principle, the a factor can 
reach ~ 1 or even 1. Before the launch of Fermi, 
several authors have argued for a strongly magnetized 
GRB central engine based on some observational evi- 
dence and its theoretical modeling. ICoburn and B oggs 
([20031 ) reported a 80% ± 20% degree of linear polar- 
ization in GRB 021206, which is best interpreted by 
synch rotron emission in a globally structured magneti c 
field ([Waxmanll2003l:ILvutikov eTaII[2003l: [Granot1[200l . 
However, a furthe r analysis of the same data cannot 
confirm the claim ([Rutledge and Fox 2004]). By model- 
ing early opt i cal fla shes within the reverse shock model, 
Zhan g et al.l (j2003) found that the early optical flashes 
with a rapidly decaying light curve (reverse shock com- 
ponent) followed by a flattening feature (forward shock) 
as observed in GRB 990123 and GRB 021211 gen- 
erally require that the reverse shock region is much 
more magnetized than the forward shock reg ion. Such 
a con clusio n was independently drawn by IFan et all 
(2002]) and[K umar and Panaitescul (|2003f ) through de- 
tailed case studies of the two GRBs. This extra mag- 
netization has to be re lated to a magnetized centra l 
engine. More recently , Kumar and McMahonl ([2008); 
iKumar and Panaitescul ([2008D analyzed the prompt 
emission data of several GRBs based on a largely model- 
independent method. They concluded that the observed 
emission cannot be produced by synchrotron emission 
within the internal shock model. They tentatively sug- 
gested a magnetized ejecta as the source of the observed 
gamma-ray radiation. 



Recent Fermi observations shed light onto the com- 
position of some GRBs. Fermi carries the Gamma-ray 
Burst Monitor (GBM) and the Large Area Telescope 
(LAT), which cover a broad spectral range of 6-7 or- 
ders of magnitude in energy. It is ideal to measure the 
broad-band spectra of GRBs during the prompt emis- 
sion phase. The first bright GRB co-detected by GBM 
and LAT, GRB 080916C, showed several nearly feature- 
less, smoothly-joint-bro ken-power- law sp ectra (the so- 
called "Band" -function, Ba nd et"al"1 (Il993[0 covering 6- 7 
decades in energy (jAbdo et al.H2009HZhang et alJ feOlO). 
Although this might not be surprising for observers (the 
Band-function has been known since the early BATSE 
years), it was somewhat surprising for modelers, since 
according to the standard fireball model, a thermal 
component associated with the fireball "photosphere" 
is expected to be very bright and should be detected 
(| Zhang and Pe 'cr 2009). This component is analogous 
to the cosmic microwave background radiation (CMB) 
associated with the hot Big Bang, and is predicted to be 
bright en ough to be detectable within a wide range of pa- 
rameters Jm eszaros and Reesjl2000tlMeszaros et al|2 002: 
iPe'er et al.l 120061: lPe'erl " l2008f ). IZhang and Pe'erl (|2009D 
analyzed this burst and argued that the non-detection of 
the thermal component strongly suggests that the ma- 
jority of energy (more than 95%) was not stored in the 
form of a "fireball" at the central engine, but was stored 
in magnetic fields which was not released until reach- 
ing a large radius. This suggests that at least for GRB 
080916C the outflow has to be Poynting flux dominated 
(PFD) with cr > (15 — 20) at the central engine an d at 
the p hotosphere as well. A follow up investigation ([Fan! 
I2010f ) on the various possibilities of hiding the thermal 
component (e.g. by invoking a smaller central engine ra- 
dius, which is not required by the minimum variability 
time scale da ta of GRB 080 916C) confirms the conclu- 
sion (see also lGao et aIll2009D 6 . All these call for a seri- 
ous re-investigation of GRB prompt emission models in 
the high-cr regime. 

In this case, magnetic energy may be sufficient to feed 
GRBs. Magnetic reconnection was suggested as a com- 
ponent for GRBs long ago (Thompson 1994). The prob- 
lem lied, however, in the intrinsic difficulty of reconnec- 
tion as it is a very slow process in ordered fields. As 
with the case for solar flares, both a slow phase of accu- 
mulation of the oppositely directed flux and a fast bursty 
phase are req uired for reconnectio n. Essential progress 
was made by (|Lazarian et al. 2003?) , who proposed a new 
scenario for GRBs by invoking self- adjusted reconnection 
based on the findin gs of fast reconnection in 3D turbu- 
lent magnetic fields ( Lazarian and Vishnia c 1999). They 
suggested that the fast bursty reconnection eventually 
occurs as a nonlinear feedback of the increased stochas- 
ticity of the magnetic field lines. The reconnection events 
start from some limited volumes and then spread in the 
form of a chain reaction as the energy is fed back to the 
turbulence and induces dramatic change in the magnetic 
field topology. The turbulent reconnection model has 

6 An alternative model is to interpret the enti re spectra as 
the emission from a dissipative fir eball photosphere ( Belo borodovl 
2010; Lazzati and BcKclman 2010). This model, although plausible 
to interpret some GRBs, is found to be difficult to interpret the 
data of GRB 080916C, see jZhang et al.ll201Ch and Sect.2.3 for a 
more detailed discussion). 
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been confirmed by recent numerical testings (Kowal et 
al. 2009). 

In this paper, we build a model of GRB prompt emis- 
sion in the high-er regime on the basis of the recent ob- 
servational and theoretical achievements. This model 
is called the Internal- Collision-induced MAgnetic Recon- 
nection and Turbulence (ICMART) model. In we 
critically review the current prompt emission models, in- 
cluding the internal shock (IS) model f £)2.1[) in the low- 
er regime, the electro-magnetic model in the extremely 
high-er regime ( ^2 .2[) and the dissipative photosphere 
model in the low- to intermediate-cr regime f£ j2.3j) . We 
argue that these models are not ideal to interpret GRB 
080916C and some other GRBs. We then delineate the 
general picture of the ICMART model in <j3l In Sj4]we dis- 
cuss the merits of the ICMART model, focusing on how 
it inherits the merits of the previous models and over- 
comes their drawbacks. In we outline several unique 
predictions of the ICMART model, which can be used to 
differentiate this model from other models based on the 
observational data. These include the existence of two 
variability components, an evolution of E p and gamma- 
ray polarization degree across each GRB pulse, as well as 
a mildly polarized external reverse shock emission from 
GRBs. We then dedicate $6] discussing how the ICMART 
model interprets the observation of GRB 080916C. The 
ICMART model is summarized in Sj7]with some discus- 
sion on its broad implications on other aspects of GRB 
physics (e.g. neutrino and cosmic ray connections with 
GRBs) as well as its possible applications to other astro- 
physical objects such as active galactic nuclei (AGNs). 
The physics invoked in this model (e.g. turbulence and 
reconnection in a high-cr, relativistic plasma) is compli- 
cated. In this paper we make the first step by delineating 
the qualitative picture of the model, and defer the more 
detailed quantitative/numerical analyses to future work. 

2. PREVIOUS GRB PROMPT EMISSION MODELS 

2.1. The internal shock model 

For a baryonic outflow (er <C 1), the standard 
GRB prom pt emission model is the internal shock 
(IS) m odel ()Rees and Meszaroslll994t iPaczvnski and Xul 
1994). It is envisioned that the GRB central engine 
launches an unsteady outflow with varying luminosity 
and Lorentz factors. Approximating the outflow as a 
series of "mini-shells" with a distribution of Lorentz fac- 
tors [r g (r min ,r max )], widths [A e (A min ,A max ) = 

(cAi m ; n , cAi max ), where At is the duration of the cen- 
tral engine activity of each minishell], and separations 
[d £ (d m in,dmax) = (c<5r. min , c<5r. max ) , where St is the du- 
ration between the end of ejecting a leading shell and 
the beginning of ejecting a trailing shell], one can ob- 
tain a series of collisions due to the interactions between 
the late, fast shells and the early, slow shells. These col- 
lisions are supersonic, resulting in internal shocks from 
which particles are accelerated and photons are released 
to power the GRB prompt emission. For two shells with 
parameters (r s ,A s ) and (Tf, Af) separated by d = cdt 
(with the fast shell lagging behind the slow shell), the 
internal shock radius is (noticing /3 = (1 — T~ 2 ) x / 2 ) 



where r.,.2.5 = r s /10 2 5 . Hereafter the convention Q s — 
Q/1Q S is adopted in cgs units throughout the text. 

2.1.1. Merits: central- engine- driven variability 

Most GRB light curves are highly variable. The in- 
ternal shock model attributes this variability to that 
of the central engine. Internal collisions are also fre- 
quently observed or inferred from other astronomical 
objects, such as pulsar wind nebula, AGN, and plane- 
tary nebulae. It is therefore natural to envision inter- 
nal collisions in GRBs as well. The connection between 
the observed GRB variability and that of the GRB cen- 
tral engine is recently strengthened by the observations 
and modeling of X-ray flares found in some Swift GRBs 
dBurrows et al.|[2005HFalcone et all 120061: iRomano et al.l 
I2006t IChincarini et al.ll2007t iFalcone et al.ll2007ft . which 
are believed to be due to late central engine activities 
([Burrows et all 120051: iZhang et all 120061: iFan and Weil 
2005 1: lLazzati and Pernal 120071 : iMaxham and ZhaneJ 
2009). A str ong support t o suc h an interpretation was 
presented by iLiang et"ail ([2006), who blindly searched 
for To of the flares based on the high-latitude "curva- 
ture" effect model of the decaying phase of the flares 
using the observed temporal and spectral data. They 
found that the required To are often near the begin- 
ning of the flares, which strongly suggests that the 
GRB central engine is restarted. Since X-ray flares and 
gamma-ray pulses share the same origin, as is demon- 
strat ed by the smooth tr ansition between the two phases 
(e.g. iKrimm et al.l [2007b . the GRB/X-ray flare data in 
general demand that the observed variability in GRBs 
should be tied to that of the central engine. The inter- 
nal shock model naturally makes such a connection, i.e. 
the observed variability time hist ory roughly traces the 
time history of the central engine (|Kobavashi et al.|[l997t 
IMaxham and Z hang 2009). Furthermore, it has been ar- 
gued that the radiation efficiency of the internal shock 
model can be much higher than that of th e external shock 
mode l to i nterpret variable light curv es (Sari and Piranl 
119971) cf. (jDermer and Mitmanl [l999h (but still not effi- 
cient enough to interpret the data, see below). This has 
made the IS model a popular theoretical model for GRB 
prompt emission for many years. 

Despite of its popularity, the IS model suffers from 
several criticisms, which we summarize in the following. 

2.1.2. The efficiency problem 

Suppose that the second, fast shell (?7i2, T2) catches up 
with the first, slow shell (mi, Ti), and that the two shells 
undergo a full inelastic collision and generate an internal 
energy U' . From energy conservation, TiTOi + T2TO2 = 
r m (mi + to 2 + C/'/c 2 ), and momentum conservation, 
TiPimi + r 2 /32W 2 = r m /3 m (mi + m 2 + U'/c 2 ), one can 
derive the Lorentz factor of the merged shell 



r — 

1 m — 



TiTOi + r 2 m 2 
toi/Ti + raijY-2 



1/2 



(3) 



The energy dissipation efficiency (which is the upper 
limit of the radiation efficiency when the "fast cooling" 
condition is satisfied) of the collision is 



Pf-P. 



2T 2 s cSt = 6 x 10 14 cm r 2 2 5 <5i_i 



(2) 



r TV 

J. 771, C 
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(Ez + Ex 



(4) 



This efficiency is typically low (jKobavashi et al.1 120071 : 
iPanaitescu et al.lll999t lKumailll999f ). The most efficient 
collisions are those with equal masses and large T ra- 
tios, but for a distribution of mass and Lorentz factor 
of a group of randomly ejected shells, the mean effi- 
ciency is low, e.g. < 15% for r max /r min < 100. For 
a standard electron equipartition parameter e e ~ 0.1, 
this model predicts an X-ray afterglow flux brighter 
than the observed flux by nea r 3 orders of magni- 
tude iMaxham and ZhansJ l2009f) . Increasing r; IS re- 
quires adjusting the Lorentz factor distribution of the 
shells, s o that the T con t rast is systematically in- 
creased (iBeloborodovl l2000t iKobavashi and SarH (20011 : 
iGuetta et al.ll200lD . There is no obvious physical reason 
why GRB central engines would satisfy these contrived 
requirements. Observationally, a detailed study of the 
GRB radiative efficiency ba sed on the Swift B AT/XRT 
data of a sample of GRBs (Z hang et al.ll2007l) suggests 
that the radiative efficiency of some GRBs can be as high 
as 90%. This is difficult to interpret within the straight- 
forward IS models. 



and the magnetic density. Comparing with the comov- 
ing dynamical time scale t' dyn = R/Tc ~ 110 s R15T2I, 
one has t' c ( , y ejP ) <C t' dyn . This suggests that electrons 
at 7 e p cool rapidly within the dynamical time, form- 
ing an electron spectrum of N(~f e ) oc 7" 2 below j e<p all 
the way to "f e . c ~ 1 <C j e ,p for the nominal parameters 
adopted. One therefore expects that the photon num- 
ber d ensity below E p follows N(E) oc E~ z l 2 ijSari et all 
1998). This is different from the typical low energy pho- 
ton spectr um observed in mos t GRBs, N(E) cx E a , with 

a 1 (jPreece et all I2000D . This is the "fast coo l- 

ing problem" of the IS model (|Ghisellini et al.|[2000T) . 
A similar argument from a diffe r ent ap proach was pre- 
sented in iKumar and McMahonl (|2008f h Since it stems 
from the synchrotron cooling argument in general, this 
problem applies to all the scenarios in which electrons 
are accelerated only once without continuous heating, 
such as the shock acceleration scenario commonly dis- 
cussed in the literature. Possibilities to alleviate this 
problem within the IS model include i ntroducing decay 
of th e shock generated magnetic fields (|Pe'er and Zhanel 
l2006f ) and introducing 2nd order stochastic Fermi accel- 
eration in the post shock reg ion due to plasma turbulence 
(jAsano and Terasaw a 2009). 



2.1.3. The fast cooling problem 

In the internal shock scenario, once an electron is ac- 
celerated at the shock front, it would cool rapidly down- 
stream through synchrotron/SSC emission. A typical 
GRB has an observed peak energy E p ~ 250 keV. This 
can be translated into an estimate of the typical electron 
Lorentz factor j ep that contributes to E p . Within the 
the standard synchrotron emission model, one has 



En 



(5) 



where V is the bulk Lorentz factor of the shocked region 
(r m in Eq. [3] for an individual collision in the IS model); 
B 1 is the magnetic field strength in the comoving frame. 
For an outflow (or "wind") with a mean kinetic lumi- 
nosity L w , the total internal energy due to shock dis- 
sipation is L w n is . Assuming that this internal energy 
is distributed to protons, electrons and magnetic fields 
in the fractions e p , e e and €3, with e p + e e + e B = 1, 
one has (L w i lls e B ) / {^R 2 cY 2 ) = B' 2 /8ir, so that TB' = 
{2L w t] 1s €b / cR 2 ) 1 / 2 . As will be evidenced soon below, the 
emission is in the "fast cooling" regime so that the en- 
ergy distributed to electrons is essentially converted to 
the observed gamma-rays. The (isotropic) gamma-ray 
luminosity (a direct observable) can be approximated as 
L 7 = L w rj is e e . Taking a typical redshift z = 1, Eq.© 
gives the constraint 



1/4 



1 + Z 



1/2 



( E p 



V 250 keV 



1/2 



(6) 



The comoving cooling time scale for an electron 
with Lorentz factor j e is t' c {^ e ,p) = {lem e c 2 ) 
/[(4/3) 7e 2 a T c( J B' 2 /8 7 rKl + y)\ = (6^m e c)/[ 7e a T B' 2 (l + 
y)] = 0.008 s l7,lB' A ~ 2 {l + y)- x , where y = U' ph /U' B 
is the ratio between the comoving photon energy density 



2.1.4. The electron number excess problem 

In the internal shock model, if all the electrons associ- 
ated with the ejecta are accelerated, for typical parame- 
ters one has too many electrons to share the dissipated 
internal energy, so that the typical synchrotron emission 
frequency is smaller than what is observed by ~ 2 orders 
of magnitude. The argument is the following. 

Since the synchrotron spectrum for the internal shock 
scenario is in the fast cooling regime, E p should be de- 
fined by synchrotron emission of the minimum energy 
electrons injected at the shock front, i.e. 7 e , p = 7e ,m. 
The mean proton Lorentz factor j p depends on the 
parameters (m, T) of the two colliding shells, and the 
relative Lorentz factor between the two shells Tf s = 
(T f /T s + T s /T f)/2. Rigorously, one can divide the shock 
interaction region into 4 regions (1: unshocked leading 
shell; 2: shocked leading shell; 3: shocked trailing shell; 
4. unshocked trailing shell). The reverse shock is typi- 
cally stronger than the forward shock, which is more rel- 
evant to GRB prompt emission. For a rough estimate, 
one has 

y p -l~ (r 43 - l)e p , (7) 

where is related to r^ s and the parameters (to, T) of 
the two colliding shells. This can be reduced to j p ~ r 43 
for e p ~ 1. With the definitions of e e and e p , one has 
7 e — 1 = (e e n p / ' e p ni){m p / m e ){^i P — 1). The minimum 
Lorentz factors of electrons and protons can be written 
as ( 7e , m -l) = <f>(p)(%-l), and ( 7p , m -l) = cf)(p)(%-l), 
respectively, with <f>{p) ~ (p — 2)/(p — 1) for p > 2, and 
cj)(p) = [ln(7A//7 m )] _1 for p = 2 (7M and j m are the 
maximum and minimum Lorentz factors of the power 
law distribution of the protons or electrons, respectively). 
Noticing HeB q /m e c = m e c 2 , where B q = 4.414 x 10 13 G, 
one can finally express Eq.([5]) into 
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; 4.4 keV 



1/2 



Hp) 



1/6 

1/2 



7-1/2 D-l 

- tJ 7,52 rl IS,14 



1 + 2 



(>-i) 2 



(8) 



Since in interna l shocks one expects ~ e 2 e < e e 
(|Medvedev|[2006D . e e < e p , and (> - 1) - 1 (Eq.0, 
unless a very large r/ s , and hence, is invoked), the 
typical E p predicted in the standard IS synchrotron emis- 
sion model is about 2 orders of magnitude lower than 
the observed value. In order to invoke synchrotron radi- 
ation as the mechanism to power GRB prompt emission, 
one needs to make an additional assumption n p /n e 3> 
1, i.e . only a small fraction of electrons are acceler- 
ated (|Bvkov fc Meszarosiri996UDaigne and Mochkovitchl 
fl998h . 

The requirement of a small fraction of accelerated elec- 
trons is also introduced in correctly calculating the syn- 
chrotron self-abs o rption frequency (v a ) in the IS model. 
IShen ancT Zhang (2009) (see their Appendix for details) 
showed that in order to get a self-consistent v a using two 
independent methods (one derived from the blackbody 
approximation in the self-absorbed regime and another 
derived from the standar d approach of applying elec - 
tron energy distribution (|Rvbicki and Lightmanl Il979l ). 
one generally requires that not all electrons associated 
with the baryonic ejecta are accelerated. 

By modeling particle acceleration in relativ istic shocks 
using particle-in-cell (PIC) simulations, ISpitkovskvl 
(2008) indeed showed that only a small fraction of elec- 
trons are accelerated, and that most electrons (~ 99% in 
number and ~ 90% in energy) are distributed in a rel- 
ativistic Maxwcllian. The numerical simulation shows 
growth of t he no n-thermal tail with simulation time. 
ISpitkovskvl (|2008l ) suspects that the thermal peak may 
eventually get "significantly eroded". In any case, it 
would be interesting to investigate the observational ev- 
idence of the putative thermal bump. The prompt GRB 
spectrum is usua lly well fit by a Band function (e.g. 
lAbdo et al.l 120091) . which does not show a thermal- like 
bump near E p (which is believed to be related to the 
injection energy of electrons due to the fast cooling ar- 
gument presented in Sect l2.1.3)) . This suggests that the 
putative thermal bump may not be significant (if it ex- 
ists at all) in the interna l shoc ks within the IS model. 
iGiannios and S pitkovskv (2009) interpret the early X- 
ray afterglow steep decay phase with significant spec- 
tral evolution (Zhang et al. 2007) as due to sweeping of 
the thermal bump across the X-ray band. This inter- 
pretation inevitably attributes the emission phase be- 
fore the steep decay phase also to the external shock 
origin. Observationally it is established that the X-ray 
steep dec ay phase is a natural extension of the prompt 
emission (iTagliaferri et al.l 120051: iBarthelmv etatl 120051 : 
lO'Brien et al.l l2006| ) . The erratic temporal behavior of 
the prompt emission is difficult to be accounted for 
within the same external shock model that interprets the 
X-ray afterglow. A more natural interpretation would be 
that the steep decay is the high-latitude emission of the 
prompt emission as long as the instantaneous spectrum 



at the end of the prompt emission phase is characterized 
by a "curved" spect rum instead of a simple power law 
(jZhang et alj|2009al) . 

2.1.5. The Amati/Yonetoku relation inconsistency 

Observationally, more energetic GRBs tend to be 
harder. This is manifested as the correlations E p oc 

E^j (jAmati et aDl2002f) an d E p oc L 7 /2 ()Wei and Gaol 
120031: lYonetoku et al.l 12004) with large scatter. The 
relations are also found to apply f or different emis- 
sion episodes within t he same burst (jLiang et al.l 120041 : 
iGhirlanda et al.l [20091 ) 7 . Inspecting Eq.©, interpreting 
this correlation within the IS synchrotron emission model 
requires that i?i4 does not vary significantly among 
bursts, i.e. the internal shock radius is insensitive to 
the GRB luminosity. Inspecting Eq.Q, this suggests 
that GRBs all share a similar Lorentz factor reg ardless 
of their luminosities (Zhan g and Meszaro s 2002a). Fur- 
thermore, increasing the average T (and hence T s ) tends 
to make a GRB softer (a larger Ru), this is in contrast 
to the naive expe ctation that high-r GRBs tend to be 
harder. Recently, iLiang et"aTI (|2010l ) discovered a tight 
correlation between T and E lj i SO based on the deceler- 
ation signature of a sample of GRBs with known red- 
shift, namely T oc E^fJ Q . Taking the trivial proportion- 
ality of L 7 oc E 1; i SO , this gives Ri$ oc T 2 oc L° 54 and 
E p oc E-°f 4 (jLiang et al.ll2010f ). which is far from the 
observed Amati relation. This is another difficulty of the 
IS model. 

2.1.6. The missing bright photosphere problem 

In order to develop strong internal shocks, the compo- 
sition of the GRB ejecta has to be baryonic, i.e. a -C 1. 
Such a baryonic outflow is believed to be accelerated from 
a fireball that initially carries most of its energy in ther- 
mal form and converts this energy to kinetic form dur- 
ing the acceleration process (lPaczvnskilll986l ; iGoodmanl 
1986; Shc mi and Piranlll990l) . As the fireball becomes 
transparent, a bright thermal emission component is ex- 
pected to leak out from the fireball photosphere, which 
forms a distinct thermal emission component in the GRB 
spectrum. Such a thermal component is predicted bright 
enough to be usually detectab le along with the IS non- 
thermal emission component ([Meszaros and Reesl 120001 : 
iMeszaros et al.l 120021: lPe'erll2008D . In the past, since the 
detector bandpass (e.g. for BATSE and Swift BAT) was 
not wide enough, there have been several speculations 
regarding this thermal component. The two leading pos- 
sibilities are: (1) the observed Band-function is the non- 
thermal component powered by internal shocks, and the 
photosphere component is either below or above the de- 
tector energy band; (2) the observed Band function is 
the superposition of a thermal component and a non - 
thermal component (|Rvdd 120051 iRvde and Pe'erll2009h . 
The excellent observational data of GRB 080916C with 
Fermi (thanks to the broad band coverage of Fermi GBM 
and LAT) show no evidence of deviation of the Band- 

7 Some arguments have been raised to show th at the global Am- 
ati relation is a pure observat i onal selection effect (Nakar an d Piranl 
l2UOirBan"d and Preec"e]|2u05t iButler et al.ll2007T) . but the fact that 
the correlation also exists internally in individual GRBs suggest a 
physical link between E p and L 7 . 
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function spectrum both below and above E p . This im- 
mediately rules out these possibilities (1) and (2). If 
one accepts that the observed non-thermal emission from 
GRB 080916C is the emission from the internal shock, a 
profound question would be "where is the thermal com- 
ponent?" . In fact, using the observed minimum vari- 
ability time scale revealed by the GRB 080916C data, 
the predicted photosphere thermal emission component 
is more than one order of magnitude brighter than the 
observed flux. This raises a severe problem to the hreball 
acceleration scenario a nd the straightforward IS model 
(|Zhang and Pe'erll2009ft . This is the "missing bright pho- 
tosphere" problem of the IS model. 

A plausible scenario is to interpret the entire Band 
spectrum as the emission from the photosphere. This 
requires energy dissipation below and above the photo- 
sphere (defined by Thomson scattering optical depth be- 
ing unity). This model is discussed in detail in Sect l2.3l 
below. 

2.1.7. Other variants of the IS model 

There are two variants of the IS model that in- 
voke different radiation mechanisms other than syn- 
chrotron emission. These variants all suffer from 
the same efficiency and fast cooling problems, but 
may introduce extra ingredients to confront other 
problems/issues discussed above. The first one 
is to invoke synchrotron self-Compton (SSC) as 
the mechanism for GRBs (IPanaitescu and Meszarosl 
20001 iKumar and Panaitescul l2008t iRacusin et alj|2008t 
Kumar and Naravanl 12009). 



length, i.e. 



space for SSC is larger than synchrotron t o interpret 
the GRB data (jKumar and McMahonll2008ft . However, 
the SSC model suffers from the following drawbacks. 
First, it generally predicts a bright prompt optical emis- 
sion com ponent. Although th is is consistent with GRB 
080319B (jRacusin et all 12008ft . prompt optical emission 
data of other GRBs are generally consistent with the ex- 
trapol ation of the gamma-ra y spectrum to the optical 
band (jShen and Zhandl2009f ). Second, a dominant SSC 
component inevitably predicts higher order SSC compo- 
nents, which greatly increases the energy budget of GRBs 
dDerishev et al.ll2001t iKobavashi et"aLl(2007tlP"iran et al.1 
|2009j). ThirdTthe synchrotron/SSC model may not re- 
produce the much more variable gamma-ray light curve 
(than the optical light curve) observed in GRB 080319B 
|Resmi fc Zhand 12010ft . Finally, E p in the SSC model 
very sensitively depends on the Lorentz factor of the elec- 
trons (E p oc 7e p ), so that a small dispersion of the ^ e:P 
distribution gives a very wid e range of E p distribution 
(|Zhang and Meszarosl l200"2ah . This makes the parame- 
ters more contrived to account for the obs erved narrow 
clust ering of E p for bright BATSE GRBs (|Preece et al.l 
12000ft . 

Another IS model variant is to invoke "jitter" emis- 
sion as the radiation mechanism (|Medvedevl I2000D . In 
this scenario, the magnetic held in the emission region 
has a random configuration, with a coherent length scale 
Xb <C Rb = jeineC 2 1 eB' , the Larmor radius. Elec- 
trons in such a field hardly make one gyration before 
the field direction changes. As a result, the typical fre- 
quency of the radiation spectrum is no longer related to 
the strength of the field, but is related to the coherent 



Ep-hT^ — (l + z)-\ 



(9) 



which is by definition much larger than E p in the syn- 
chrotron emission model (Eq.[5]) given the same Te.p- 
This eases the constraint discussed in 32.1.41 However, 
since there is no prediction on As from the first prin- 
ciple (and it is not clear whether such a characteris- 
tic scale exists), it is not easy to assess how E p de- 
pends on model parameters. In particular, the corre- 

1/2 

lation E p oc L 7 established in the synchrotron model 
(Eq.[Sl ) is no longer straightforw ard to justify. Re- 
cently, [Shonrand Spitk ovskvl (|2009t ) synthesized the par- 
ticle spectrum from PIC simulations and concluded that 
the spectrum is entirely consistent with synchrotron ra- 
diation in the magnetic fields generated by Weibel insta- 
bility The "jitter" regime is recovered only when one 
artificially reduces the strength of the electromagnetic 
fields. 

2.2. The electromagnetic model 

In another extreme, i.e. a ^> 1, a so-called Electro- 
Magnetic model (EM model) (|Lvutikov and Blandfordl 
120031 : lLvutikovll2006bft has been proposed. In the high-c 
regime, the comoving Alfven s peed is close to speed of 
light, which can be written as (|Jacksonlll975ft 



(c 2 + n.NR 2 ) 1/2 ' 



The allowed parameter where 



^A,NR 



B' 



V47T/3' 



(10) 



(11) 



is the comoving Alfven speed in the non-relativistic 
regime (noticing the definition of a in Eq.pQ), one can 
write the Lorentz factor of the comoving Alfven wave as 



7'a 



(1 + *) 



1/2 



(12) 



For a cold plasma (i.e. the comoving thermal energy 
density much smaller than the rest mass energy den- 
sity, and hence, the magnetic energy density), gas sound 
speed is much less than the Alfven speed. Therefore 
"f' A is also approximately the comoving Lorentz fac- 
tor of other magnetoacoustic waves. The EM model 
(|Lvutikov and Blandford 2003) applies to the "sub- 
Alfvenic" regime, i.e. T < Y A = (1 + c) 1 / 2 , or roughly 

(13) 



a > a, 



= T 2 - 1 



io 6 r2 



at the deceleration radius. According to this model, 
the Poynting flux is dissipated through electromagnetic 
current-driven instabilities. The model has some novel 
features that are not shared by the fireball IS model. 
For example, the emission radius where strong magnetic 
dissipation occurs is R ~ 10 16 cm, much larger than 
that of internal shocks defined by the minimum vari- 
ability time scales (Eq.[2]). This large radius is consis- 
tent with independent co nstraints on GRB emission site 
using different methods (iLvutikovl l2006at IKumar et all 
20071: iShen and Zhand 120091: iGupta and Zhand 120081 : 
Zhang and Pe'erl 120091 ) . The EM model also justifies a 
"structured" jet with energy per unit solid angle drop- 
ping with angle as E(6) oc 6*~ 2 , as has been invoked 
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in some GRB phenomenological models jM eszaro s et al.l 
1998; iRossi et alJ[2002l IZhang and Meszarosll2002b[ ) . 

There are however two major issues related to the EM 
model. The first one is regard ing its extremely high-cr 
value at the deceleration radius. iLvutikov and Blandford 
(2003) argued a force-free flow by invoking a ~ 10 9 
beyond the photosphere. However, considering more 
realistic models to launch a magnetized outflow from 
the central engine wou ld lead to a rang e of a values 
that are below a c (e.g. ISpruit et al.ll200l1) . At the ex- 
treme, even the outflow consists of cold (pressureless) 
matter accelerated exclusively b y the magnet ic field, 
the achievable a is at most a c (jMichell I1969I ). Such 
a maximally magnetized outflow is achievable only at 
a specified geometry, i.e. a purely radially expand- 
ing outflow. For more general geometric configurations, 
the achievable a should be below a c , and could be 
as lo w as a < 1 (Beg elman and Lil Il994t ISpruit et al.l 
12001 . If one consider s the "collapsar" scenario of lon g 
GRBs (|Wooslevl [1991 iMacFadven and Wooslevl fl999h . 
a PFD jet launched from the central engine needs to 
penetrate through the stellar envelope. The a value 
would be further degraded. Recent numerical simula- 
tions (jTchekhovskov et al.|[200l lNagatakill2009ft suggest 
that a high-cr flow at the base of the flow would become a 
mode rate-cr flow as it escapes the star (cf. IWheeler et al.1 

Another issue of the EM model is related to GRB 
variability time scales. According to the EM model, 
the GRB central engine is assumed to launch a mag- 
netic bubble, which expands and exits the star with- 
out degrading a, and magnetic energy is finally dissi- 
pated by the c urrent-driven instability at a la rge radius 
R ~ 10 16 cm (jLvutikov and Blandford! I2003D . The ob- 
served rapid variability of the GRB light curves is in- 
terpreted as due to emission from some "fundamental 
emitters" that are moving relativistically with respect to 
the dissipation region, which itself is moving relativis- 
tically with respect to the observer. A related model 
invoking relativistic, possibly magnetic tur b ulenc e was 
recently proposed by iNaravan and Kumarl ()2009l ) and 
iKumar and Naravanl (|2009[ ). Within such a picture, the 
observed variability time scale does not reflect the be- 
havior of the GRB central engine. As discussed earlier in 
£12.11 the connection between at least the slow variability 
time scales with GRB c entral engine activ ity is strongly 
supported by the data ([Liang et al.ll2006l ). especially in 
view of the discovery of X-ray fla res that extend cen- 
tral engine activity to later times dBurrows et al.l [20051 : 
IChincarini et al.l 120071 : iFalcone et al.ll2007f ). Attributing 
the observed variability to random motion of the rela- 
tivistic "fundamental emitters " in the emission region 
(jLvutikov and Blandford 2003]) would then require argu- 
ing th at the "Xh-reset" e ffect of the last pulse and X-ray 
flares ([Liang et al.ll2006l ) is purely due to a chance coin- 
cidence. 



2.3. MHD and dissipative photosphere models 

If the magnetized ejecta can be approximated as a 
fluid and described by magnetohydrodynamics (MHD) 
(|Spruit et al.ll200lD . the GRB models are in the MHD 



regime 8 . Within this regime, outflows launched from the 
central engine would carry a dynamically important mag- 
netic field, but it can be still approximated as a MHD 
fluid 9 . The a value may decrease with increasing radius 
due to two mechanisms. First, below the photosphere, 
reconnection dissipation energy is thermalized in the 
plasma because of the high optical depth of the radiated 
photons. The thermal energy is then converted to bulk 
kinetic energy sim ilar to the fireball acceleration process 
(|Drenkhahnl 120021: iDrenkhahn and Spruitl 12002ft . This 
scenario is more relevant to the so-called "striped wind" 
magnetic field configuration with alternating magnetic 
polarity (relevant to pulsar-like central engine with mis- 
aligned magnetic and rotational axes). For helical mag- 
netic configurations (relevant to black hole central engine 
with rotational and magnetic axes aligned), the neighbor- 
ing ordered field lines typically have the same orientation, 
so that reconnection is greatly suppressed. Second, the 
Poynting flux may be direct ly converted to kinetic en- 
ergy without dissipation (e.g. IVlahakis~an d Konigli [20031 : 
iKomissarov et all 12009) under the pressure gradient of 
magnetic fields in an MHD flow. The efficiency of such 
a conversion is, however, subject to uncertainties. With- 
out an external pressure confinement, the Poynting flux 

energy cannot be fully converted to kinetic energy, and 

1/3 

the flow can be only accelerated up to r + „ t , where r t , n t 
is th e Lorentz factor for total conversion ([Beskin et all 
1998). This was the origin of the well known "<r prob- 
lem" of the Crab nebula. The jet acceleration can be ef- 
ficient with an external confinement , which is relevant to 
GRBs in the collapsar scenario (e.g. lTchekhovskov et al.l 
2008). In any case, the pressure gradient has to be small 
enou gh for efficient acceleration to happen ([Lyubarskyl 
2010). Overall, if the jet has ^ > 1 at the central en- 
gine, the a value at the GRB emission region can range 
from a < 1 to a < o-q, depending upon how efficient 
the conversion from P oynting flux to kinetic flux would 
be. iLvubarskvl (|2010l ) showed that it is difficult to have 
a completely matter dominated jet (say, a < 0.1) at the 
GRB emission radius for a dissipation- less jet with an ini- 
tial high <7, even if the most efficient conversion occurs 
(see also Levinson 2010). 

8 ILvutikov! 11200611) defines this regime as 1 < a < cr c = 
r 2 — 1 ~ 10 6 r|, opposed to the sub-Alfvenic regime invoked by 
the EM model. He also relates the "super-Alfvenic" and "sub- 
Alfvcnic" regimes to whether or not a shock can exist in the 
ejecta. Such a connection is relevant only for ID ejecta-medium 
interaction. A more relevant condition to discuss shock forma- 
tion condition in the high-cr regime is whether the magnetized 
eject a encounters a stronger pressure than its own ma gnetic pres- 
sure (Zhang and Kobayashi 2005; Mizuno et al. 2009a). For colli- 
sions between two shells with the same cr value, internal shocks can 
develop even for a small relative Lorentz factor between the two 
shells. For an expanding shell with small width (thin shell) and 
conical geometry, the reverse shock cannot form for a > (10 — 100) 
(Zhang and Kobayashi 2005). 

9 Strictly speaking, the outflow may not satisfy the MHD con- 
dition at all radii. Since the particle density n oc R~ 2 decays 
with radius faster than the Goldreich- Julian density n GJ oc R , 
the MHD condition would break above a critical radius R MHD . 
If this radius is be yond the deceleration radius (R MHD > ^?dcc) 
(ISpruit et al.l 1200 J l. then the MHD condition would not break 
before afterglow sets in. However, if R MHD < Rdeo the MHD 
condition breaks at R MHD and a strong magnetic dissipation oc- 
curs between R^„^ and -Rhrh- This happ ens if a is larger than a 
few hundred (Zhane and Mcszaros 2002a). We do not discuss this 
regime in this paper by limiting our discussion to the regime of 
1 < a < 100. 
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In the past, th e GRB prompt emission models in 
the MHD regime (ISpruit et ahl 120011: iDrenkhahnl [2001 
Drenkhahn and Spruitll2002t IVlahakis and Konig]||200a 



Giannios 2008) have been "semi-fireball" like. It is gen- 
erally assumed that an initial high a is reduced to a < 1 
before the jet reaches the internal shock radius, so that 
the traditional IS model can still apply at i?is- Such a 
model shares essentially the same advantages and draw- 
backs of the IS model. The data of GRB 080916C give 
a strong constraint on such a model: the lack of a pho- 
tosphere component requires that either a 3> 1 at the 
photosphere, or a < 1 already at the photosphere. For 
the former, the magnetic-to-kinetic energy transition ef- 
ficiency is low. One cannot have matter-dominated in- 
ternal shocks. For the latter, the magnetic acceleration 
proceeds in a "cold" and efficient way, which is difficult 
to achieve theoretically. The external confinement of the 
high-er flow (which is required for efficient magnetic ac- 
celeration) would inevitably cause heating of the jet to 
result in a bright photosphere, which is constrained by 
the data of GRB 080916C. Even if the flow could be- 
come baryon-dominated through magnetic acceleration 
without a hot photosphere component, the internal shock 
efficiency is so low that only a small fraction of energy 
(say ~ 10%) is converted to gamma-ray radiation. This 
greatly increases the total energy budget of the outflow, 
and demands an even higher a at the central engine to 
begin with (from the missing bright photosphere argu- 
ment). This makes it even more difficult to efficiently 
convert Poynting flux energy to kinetic energy. 

A plausible scenario is to argue that the observed 
spectrum is dominated by the emission from the pho- 
tosphere itself, and that the internal shock emission 
is too weak to be detected. Within such a sce- 
nario, there is no "missing photosphere" problem, 
since the observed spectrum is the photosphere emis- 
sion itself. One then needs to argue that the en- 
tire non-thermal GRB spectrum is the reprocessed pho- 
tosphere thermal emission. This requires significant 
energy diss i patio n below and above the p hotosphere 
([Thompson! 119941 iRees and Mes zaros 20 051: iThompsonl 
' 20061 iThompson et all l2007t IGhisellini et al.l 120071: 
Gianniosl2008HBeloborodovl2010l : lLazzati and Begelmanl 



20101 ). Within the MHD models, the dissipation source is 
continuous magnetic reconnection (for the striped-wind 
geometry). Such a model can be also developed for a 
neutron-rich baryonic flow with a <C 1, in which neu- 
tron dec ay provides a sour ce of continuous heating in 
the jet (jBeloborodovi [2010). By properly considering 
Compton upscattering of photosphere thermal photons, a 
hard power law tail can be p roduced (Bclobo rodovll2010l : 
lLazzati an d Begclman 2010). This model can also nat- 
urally interprets GRB variability in terms of the central 
engine energy injectio n history. The radia tive efficiency 
is also naturally high (jLazzati et al.ll2009T ). Confronting 
the GRB 080916C data, however, this model faces two 
difficulties. 

First, the maximum photon energ y in the dissipative 
photosphere model is ~ 0.1 — 1 GeV ()Beloborodovll2010l ; 
Lazzati and Begclman 2010). The detection of 13.2 GeV 
photon associated with one of the GBM pulses in GRB 
080916C (which corresponds to 70.6 GeV rest-frame) 
then disfavors the photosphere origin of gamma-ray emis- 
sion in GRB 080916C. It has been argued that the GeV 



emission of LAT GRBs inclu ding GRB 080916C may be 
of th e external shock origin (jKumar and Barniol Duranl 
120091: IGhisellini etall 12010ft . However, the required 
forward shock parameters are extreme, which are 
not easy to accomm odate within the known rela- 
tiv ists shock models ([Kumar and Barniol Duranl 120101 : 
|2010bl iPiran and Nakarl l2010ft~ Furthermore, the 
LAT-band emission and GBM-band emission generally 
trace each other during the prompt emission phase 
(jZhang et al.l [2010h 10 . For GRB 080916C, the time- 
resolved spectral analysis with the finest temporal resolu- 
tion defined by statistics suggests that the Band function 
fits well the data for every time bin throughout the burst. 
The peak of the GeV emission in the logarithmic light 
curve coincides with the second peak of the GBM light 
curve. All these strongly suggest that the entire emission 
of GRB 080916C is from the same emission region, and 
is therefore likely of an internal origin (see Zhang et al. 
2010 for more detailed discussion). The long term GeV 
emission, on the other hand, decays slower than the MeV 
emission. It may be dominated by a different emission 
component (e.g. external shock). 

Second, although the high energy photon spectral 
index above E p can be well reproduced in the dis- 
sipative photospher e model by Compton upscat tering 
(jBeloborodovi 120101 : lLazzati and Begelmanl [2010), the 
low energy photon spectral index below E p is typi- 
cally much harder than what is observed. For up- 
Comptonization of a blackbody spectrum (e.g. within 
the context of soft gamma-ray repeaters and anoma- 
lous X-ray pulsars), the Ray leigh- Jeans low energy spec- 
tral in dex (Band-function index a — +1) is hardly af- 
fected (jNobili et al.ll2008T )). By invoking dissipation be- 
low the photosphere (e.g. neu tron heating) , this index 
can be modified to a — +0.4 ( Belo borodovl [2010l) . still 
much harder than the observed typi cal value a ~ — 1 
(|Abdo et all 120091 : IZhang et all 12010ft . One possibility 
would be to assume that the observed low energy spec- 
trum is the supe rposition of the photosp here emission of 
many shells (e.g. lToma et al.ll2009L 12010 ) . Although this 
is not impossible, it requires properly arranging the lumi- 
nosity and Lorentz factor of many mini-shells to mimic 
a typical a = — 1 Band spectrum. A detailed time- 
resolved spectral analysis of GRB 080916C and other 
"Band-only" GRBs suggests that the low energy photon 
index a remains essentially u nchanged as the t ime bin 
becomes progressively small (jZhang et alJl20foh . This 
at least disfavors the possibility that the observed time- 
integrated spectrum is the temporal superposition of the 
photosphere emission of many shells. One plausible solu- 
tion would be to introduce synchrotron and synchrotron 
self-Compton (SSC) of the electrons. In the MHD model, 
if continuous magnetic heating is operating, the SSC 
component would produce a typical Band spectrum in 
the MeV range, if most energy dissipation happens in the 
T homson-thin reg ion above the photosphere (e.g. Fig. 2 
of iGianniosI 12008 ) . Such a model predicts a bright opti- 
cal emission component (synchrotron). A lthough this is 
consi stent with the case of GRB 080319B (|Racusin et al.l 
2008) , it is inconsistent with most of the prompt optical 
emission data or upper limits, which are consistent with 

10 This is the case even for GRB 090902B, which clearly shows 
two distinct spectral components. 
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or below the ext rapolation of gamma-r ay spectrum to 
the optical band (jShen and Zhandl2009| ). 

2.4. Summary 

In summary, the new Fermi LAT data of GRB 080916C 
raises a challenge to the traditional fireball IS model, 
which has some advantages and drawbacks. The EM 
model invokes too high a a at the deceleration radius, 
which may not be achieved in nature. We consider the 
modcratc-cr MHD regime a relevant regime for GRB out- 
flows. However, the current semi-fireball MHD models or 
the dissipative photosphere models are not ideal to ad- 
equately interpret GRB 0809I6C. The ICMART model 
proposed in this paper is meant to inherit the merits of 
the existing models and to overcome their drawbacks. 

Since most LAT GRBs are found similar to GRB 
080916C to show "Band-only" time-resolved spectra 
(jZhang et all 12010ft . we argue that the proposed model 
would be relevant to most GRBs. 

3. THE ICMART MODEL 

3.1. Basic assumptions 

The Internal-Collision-induced MAgnetic Reconnec- 
tion and Turbulence (ICMART) model is based on the 
following two assumptions: 

1. The GRB central engine is intermittent in nature, 
which ejects an unsteady outflow with variable 
Lorentz factors and luminosities but a nearly con- 
stant degree of magnetization. Approximating the 
unsteady outflow as some discrete shells with vari- 
able Lorentz factors, these magnetized shells collide 
with each other at the conventional internal shock 
radius Ri S (Eq.0). 

2. The GRB central engine is highly magnetized, 
which ejects a high-er flow from the base of the 
central engine. Although various mechanisms (e.g. 
magnetic acceleration and reconnection) may re- 
duce a as the outflow streams outwards, the ejecta 
is still moderately magnetized in the GRB emission 
region, with a ranging from 1 < a < 100. 

In an ideal magnetically driven MHD jet, magnetic en- 
ergy is converted into kinetic energy during the expan- 
sion. Usually it is assumed that do > 1 at the cen- 
tral engine with Tq ~ 1. The magnetization param- 
eter a t hen decreases with radius as F increase s with 
radius dBeskin et all Il998t IVlahakis and Konigll 120031: 
Komissarov et al.l 120091: iLvubarskvl 120101 : iGranot et al.l 
20101) . Without external pressure confinement, the early 

acceleration phase reaches a Lorentz factor F ~ o~^ 3 ', 

2/3 

with remaining magnetization a ~ a . Beyond this 
radius acceleration is slow and weakly dependent on R, 
e.g. T oc i? 1 / 3 , and a cx iT 1 / 3 (|Granot et al.l[2010h . If 
(To remains roughly constant, one expects a nearly con- 
stant cr at large radii, but in the meantime does not 
expect a large variation in L. We note, however, that 
near the rapidly-rotating magnetized central engine, ef- 
ficient acceleration of ejecta may occur below the light 
cylinder thanks to direct electric field accelerate mech- 
anisms (e.g. frame-dragging-induced deviation from the 
force-free condition, Muslimov & Tsygan 1992 within the 



context of pulsars). As a result, at the light cylinder (ef- 
fectively the base of central engine), the ejecta may have 
Lo ^> 1. The rapid variation of the jet power (e.g. due 
to variable accretion rate or variable spindown rate of 
the central object) may lead to a variable Fo and some- 
what variable <7o at the central engine. Considering later 
magnetic acceleration, one may get a large T variation 
within the outflow, with er not significantly varying at a 
same radius. Our above two assumptions may be then 
fulfilled. 

3.2. Basic parameters of the GRB ejecta plasma 

In preparation of proposing the model, it is informative 
to summarize the basic parameters of GRB ejecta. The 
gamma-ray emission radius R can be in principle in the 
range of 10 11 — 10 17 cm (bracketed by the photosphere 
radius and the external shock deceleration radius). As 
will become evident later, several collisions are needed 
to trigger an ICMART event, so that one would have 
-Ricmart > Ris- In the following, in order to make dif- 
ferentiation with the IS model, we will normalize all radii 
relevant for the ICMART model to R = 10 15 cm R 15 , and 
denote R as R 1CMART whenever relevant. 

1. Length scales: At i? ICMART , the "thickness" of the 
ejecta is related to the radius R 11 , i.e. 

A' = — ~ 3.2 x 10 12 cmR 15 T^l (14) 

in the comoving frame, and 

A = j| ~ 10 10 cimR^ 2 , (15) 

in the rest frame of the central engine (lab frame). 
Assuming a conical jet with the opening angle 0j , 
the cross section radius of the emission region is 

Re = R9 3 = 8.7 x 10 13 f |m R 15 cm (16) 

in both the lab frame and the comoving frame. 
Typically one has R s > A' > A. GRB ejecta 
are therefore also considered as "flying pancakes" . 

2. Plasma number density: In a conical jet, density 
drops with radius as R~ 2 . For a hydrogen ejecta 
with a total "wind" luminosity L w , Lorentz factor 
r, and magnetization parameter er, the comoving 
ejecta proton number density is 

n == 

p 4n(l + o-)R 2 r 2 (m p + Ym e )c 3 

~ 1.8 x 10 7 cm' 3 L^T^R^m- 1 ^ + a)^ 1 (17) 

where m = 1 + Ym e /m p ~ 1 if the lepton (pair) 
multiplicity parameter Y <C m p /m e . In the rest 
frame of the central engine, the ejecta proton num- 
ber density is 

n p = Tn' p ~ 5.6 x 10 9 cm" 3 ^ ^iR^m^il + a)^ 1 . 

(18) 

11 This approximation is based on the assumption of shell- 
spreading. FOThj i gh ; o_sheJls, the spreading condition may be more 
stringent (Granot et al. 2010). The shell width is then related to 
the duration of the central engine, not a function of R. In any case, 
the following estimates are valid to order of magnitude. 
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The electron number (or more generally the lepton 
number) densities are 

n' e — Yn' p , n e — Yn p . (19) 

3. Magnetic field strength: The comoving magnetic 
held strength is 

/ T \ V 2 

/ L w r; \ 



B' 



\T 2 R 2 cl + o 



1/2 



-1.8 x 10 3 G L^T-iR-i . (20) 

In the lab frame, the magnetic field strength is 

B = TB' ~ 5.8 x 10 5 G (j^j ' L% 2 R^ . (21) 

This B field is accompanied by an induced E = 
-V x B field for an ideal MHD fluid. 

4. Collisional mean free path and time scale: For 
Coulomb collisions, the strong collision radius may 
be defined by e 2 /r co i — kT so that r co i ~ e 2 /kT ~ 
(1.7x 10~ 3 cm)/T. Here kT denotes more generally 
the typical energy of the particles (which are not 
necessarily in thermal equilibrium) . The comoving 
collision mean free path of electrons can then be 
estimated as 

l'e,coi = {<*r 2 col )- 1 ~6.5x 10 17 cm 

x L-)^l.AmY-\l + a) x Tl w . (22) 

In order to have the plasma in the "colli- 
sional" regime, one needs to require l' e col < 
A'. This is translated into T e < T e c = 

2.2x 10? K L^T-rR^m-^Yy^l+ar, 112 . 
Downstream of the GRB internal shocks, the effec- 
tive proton temperature is of the order of T p — 
(T ud - l)m p c 2 /k ~ 1.1 x 10 13 K {T ud - 1) (where 
T uc i is the relative Lorentz factor between upstream 
and downstream), and the electron temperature is 
even higher by a factor of (e e /e p )(m p /m e ). As 
a result, GRB shocks must be "collisionless" . In 
the upstream, the temperature can be lower, but 
since the flow is relativistic, any inhomogeneity 
in the velocity field would result in "heating" in 
the flow. Other processes (magnetic reconnection, 
neutron decay, etc) would also enhance heating, 
so that the electron temperature may be main- 
tained to be significantly above T e c . For a rea- 
sonable estimate, the electrons in a relativistic 
flow would have at least a relativistic tempera- 
ture T e — m e c 2 /k = 5.9 x 10 9 K. The collision 
time scale can be estimated as t' co1 — l' ecol /v' e . 
For a non-relativistic temperature, one has v' e = 
(2kT e /m e ) 1 ^ 2 , so that the comoving collision time 



'col, NR. 



_ 1/2 



TOe /2 (fcr e ) 3 / 2 /(%/27re 4 <). In the GRB 
ejecta (even without shock heating), it is very pos- 
sible that the electrons have a comoving relativistic 
temperature, so that v' e ~ c. The comoving colli- 
sional time is therefore 



'col,R - 



e,co! 



2.2 x 10 7 s 



x K^l.AmY-^l + a^T 2 , (23) 



which is > t' dyn = R/Tc - 110 sR 15 T^l, the 
comoving dynamical time scale. This again sug- 
gests that ejecta is collisionless even without strong 
shock heating. 

5. Gyroradii and gyrofrequencies: Another relevant 
length scale is the particle gyro-radius in magnetic 
fields. Without direct collisions, the GRB ejecta 
can be still approximately described as a "fluid" 
macroscopically. This is because particles are inter- 
acting with each other through magnetic fields mi- 
croscopically. Fundamentally the smallest length 
scale is defined by particle gyration. The comoving 
frame gyro(cyclotron)-radii are 

1/2 



' B,e 



B , - 0-93 cm j e L w J, 2 r 2 . 5 i?i 5 



(24) 

for electrons (where 7 e is the electron Lorentz fac- 
tor), and 



B,p 



~f P m p c 
eB' 



1.7xl0 3 cm 7p L tUi 1 5 / 2 2 r 2 . 5 -Ri5 / ' ' " 



1/2 



(25) 

for protons (where j p is the proton Lorentz fac- 
tor). Given the typical j e values to interpret the 
GRB emission (Eq.[B]) and the related 7 P (which is 
typically smaller by a factor of e p m e / e e m p ), both 
radii are < A'. If the typical viscous length scale is 
not much larger than the gyration radius, the fluid 
description of the GRB plasma is justified. The 
comoving gyrofrequencies are 

1/2 



J B,e 



eB' 

m e c 



3.2 xio 10 B-XM^Mi^: 



for electrons, and 
eB' 



(26) 



J B.p — — w B,c 

m v c m v 



.1.7 x 10 7 s-'L^T^R^ 1 



1 + cr 



1/2 



(27) 



for protons. Both are much larger than the inverse 
of the comoving dynamical time, i.e. (R/Tc)^ 1 ~ 
9.5 x 10~ 3 s _1 r 2 .5i?j~ 5 1 , again justifying the fluid 
description. 

6. Plasma frequencies and plasma skin depths: The 
comoving relativistic plasma frequencies are 

for electrons, and 



4irn p e 
%m p 

-1/2 7-1/2 



/„2\ 



5.5 x 10 6 s" 1 



x % 1,2 L%T- 2 \R^m-^{l + CT )r 1/2 (29) 

for protons. Here 7e and 7 p denote the mean 
Lorentz factor of the relativistic electron and pro- 
ton gas, respectively. The corresponding plasma 
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skin depths are 
8[ - 130 cm,/ 2 Y 1 2 L 



,1/2/ 



52 1 2.5-Kl5»W"(l+0-)i 
(30) 



1/2 



for electrons, and 
= * 5.4xl0 3 cm7i/ 2 L; 1 /2 2 r 2 .5i?i5m 1 / 2 (l+ ( 7) 1 /2 

(31) 

for protons. Physically these quantities are rele- 
vant to a weakly magnetized ejecta with random 
B fields (u <C 1). For a cr > 1 flow with ordered 
magnetic fields, the plasma oscillation frequencies 
and skin depths are relevant only in the direction 
parallel to the magnetic field lines. In the perpen- 
dicular direction, the gyrofrequencies and gyroradii 
are more relevant. 

3.3. Turbulent nature of the GRB ejecta 

Turbulence is believed to be ubiquitous in astrophysi- 
cal systems. This is because the Reynold's number, the 
ratio between the viscous diffusion time t u = L 2 jv and 
the relative flow time scale ry = L/5V, 

R e =^, (32) 
v 

is ^> 1 (mostly because of the large L involved in astro- 
physical systems) , where L and SV are the characteristic 
length and relative velocity of the flow, and 



v ~ c s l 



(33) 



is the kinematic viscosity as defined in the MHD equation 
of motion 



,M7sr + v.w) 



-VP + J x B + pi/V 2 V , (34) 



c s is sound speed, and I is the mean free path of micro- 
scopic interactions that define the viscosity. 

Since shear motion tends to distort the fluid while the 
viscous term tends to smear the distortion of the fluid, 
the flow would become highly distorted and turbulent 
when the Reynold's number is 1. 

Being magnetically dominated, most kinetic motions 
are concentrated in the direction perpendicular to the 
magnetic field. The corresponding perpendicular viscos- 
ity is (Spitzer 1962) 



1.7 x 10 _2 cm 2 s" 1 n 



In A/ (VTb 2 



(35) 



where In A is the Coulomb logarithm. A = 
3/2eV/c 3 T 3 /7mmin(l, ^4.2 x 10 5 /T). For typical pa- 
rameters adopted in this paper, one has v± ~ 3.7 x 
lCT 4 cm 2 s- 1 < L5V ~ 10 24 cm 2 s" 1 , so that R e ~ 
10 28 3> 1. Therefore th e GRB ejecta is turb ulent in na- 
ture. In the low-cr case, IZhang et all (|2009b[ ) has shown 
numerically that the GRB ejecta with a mildly relativis- 
tic relative motion (relevant to internal shocks and late 
external shock in the trans-relativistic regime) quickly 
turns turbulent if a Kelvin-Helmholtz instability is trig- 
gered. In the high- cr regime, due to the strong magnetic 
pressure in the ordered magnetic fields, the condition to 
trigger turbulence would be more demanding (see §3.5j) . 
In any case, strong anisotropic turbulence can develop 
once the triggering condition is satisfied. 



For a resistive magnetized flow, another relevant di- 
mensionless parameter is the magnetic Reynold's num- 
ber, the ratio between the magnetic resistive diffusion 
time Tdif = L 2 /rj and the flow time Tt = L/SV, 



LSV 

1 



(36) 



where rj is the magnetic diffusion coefficient, which is 
defined in the diffusive MHD induction equation 



— = V x (VxB)+i)V 2 B 

at 



(37) 



For the reason described above, we adopt the perpendic- 
ular resistivity in a strong magnetic field (Spitzer 1962), 



1.3 x lCFcmV 1 



ZlnA 

T 3 / 2 



(38) 



For the typical parameters adopted in this paper, one 
has r]± ~ 1 cm 2 s _1 . This gives a huge R m number, 
R m ~ 10 24 . The maximum resistivity is the "Bohm" 
diffusion, i.e. 

m < r B V ~ r' B e c. (39) 



This gives 



R 



m.B 



A>/r> B ~3AxlO^L^ 2 



1/2 



(40) 

This is still a large value, suggesting that magnetic fields 
can be also highly distorted and turbulent if the turbu- 
lence triggering condition is satisfied. 

Unlike the hydrodynamical turbulence that displays a 
Komolgorov solution, i.e. dE{k) oc k~ 5 / 3 (where k is the 
wave number and E(k) is the energy per unit wave num- 
ber), MHD turbulence is anisotropic and has different 
scaling in the directions perpendicular to and along the 

field line, namely E[k±) oc k ± (Komolgorov-type) and 
E{k\\) oc k^ 2 (|Goldreich and Sridharl fl995l iCho et all 

2002). While the kinetic power in the turbulence drops 
with k quickly, the power of magnetic fields does not 
drop significantly with k. As a result, the "eddies" in 
smaller scales are even more stretched and appear elon- 
gated along the local magnetic field. The MHD turbu- 
lence in the a > 1 regime is not well studied, but qualita- 
tively it should be the natural e xtension of the ph ysics in 
the low /3 = Pgas/Pmag regime (|Cho et al.ll200l 12 . The 
main extension is that the relative motion must be rela- 
tivistic in order to distort the magnetic field lines. This 
is not very demanding in a GRB, since the outflow itself 
is highly relativistic. The MHD turbulence in a GRB 
ejecta is relativistic in nature. 

3.4. Magnetic reconnection in the GRB ejecta 

Reconnection plays a fundamental role in both labora- 
tory and astrophysical settings, although the full details 
of the process was not well understood until recently. 
The main difficulty was that in a stea dy state, the recon- 
nection process proceeds very slowly ( Sweet 1958; Parker 

12 Notice that /3 is not simply since P ga s is the thermal 

pressure of the gas, while in the definition of a (Kq.[T]), the rest 
mass energy density is invoked. Even in the very low-/3 regime 
(Pgas ~ 0), the a value can be still below unity. 
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1957), and is not adequate to account for the abrupt 
reconnection events observed in the lab and various as- 
trophysical phenomena, such as solar flares. This was 
merely a theoretical problem as in reality reconnection is 
indeed fast since astrophysical magnetic fields would be 
completely entangled due to turbulent motions. In the 
standard non-relativistic Sweet-Parker scenario, two sets 
of field lines with opposite orientations approach each 
other and reconnect within a layer of thickness S and 
length C, which satisfy 



8 _ V in _ c-1/2 



where the Lundquist number 



S 



CV A 
V 



(41) 



(42) 



is essentially the magnetic Reynold's number in the 
Bohm diffusion approximation, with SV replaced by the 
Alfven speed Va, and L replaced by C. Notice that 
for individual reconnection events, it can well be that 
C <C L ~ Rg (Eq.[IB]). In any case, S is usually a 
large number, so that the Sweet-Parker process is an 
extremely slow process, i.e. Vi n = VaS~ x I 2 <C Va- 
A fast steady-state reconnection scenario was proposed 
bv iPetschekl (|196 4). which invokes a shorter width C of 
the resistive layer. However, such a scenario is unsta- 
ble unless the magnetic diffusion r\ increases near the 
X-point. Simulations with a constant r\ suggests that 
an initial Petschek configuration would quickly revert to 
the S w eet-Parker configuration (lUzdenskv and Kulsrudl 
l2000t ). lLazarian and Vishniad (|1999| ) proposed that re- 
connections in magnetic fields with stochasticity can pro- 
ceed rapidly, thanks to the turbulent nature of the mag- 
netized fluid that both broadens the reconnection zone 
and allows many independent flux reconnection events to 
occur simultaneously. This turbulence model of 3D re- 
connection is confirmed by numerical simulations (Kowal 
et al. 2009), and is the basis of our discussion of the IC- 
MART model 13 . 

In the presence of turbulence, a magnetic field recon- 
nects over a local scale An rather than the global scale 
C. This is the scale over which a magnetic field wanders 
away from its original direction by the thickness of the 
Ohmic diffusion layer. Accordingly, it is the parameter 



MVa 



(43) 



rather than the Lundquist number S that matters. The 
comoving local reconnection rate is given by V^ ec loc w 

V A s~ 1/2 < Va ~ c (|Lazarian and Vishniad 1 1 9991 ). 

One may estimate the required global reconnection 
speed to power a GRB. For the high-cr flow invoked 
in this paper, a significant fraction of the magnetic 
field energy is converted into radiation. This re- 
quires that the dissipated magnetic energy in a shell 
with comoving width A' within the time scale At' 
can power the observed gamma-ray luminosity, i.e. 
T 2 (B' 2 /8n)4nR 2 A'/At' ~ 10 52 £ 7 , 52 . To compare 

13 We note that the study of reconnection i s an rapidly evolving 
field IIDrake et al,ll200H : ILoureiro et al . 2007, 2009; Z enitani et"afi 



against the definition of B' (Eg .[20]). this demands that 
the global reconnection speed 



^rec, global * ■» 



At 7 



£7,52 1 + cr 
L w ,52 a 



(44) 



In turbulent reconnection, the global reconnection rate 
is boosted from VJ.' CC loc by a factor of ~ L/\\\, 
since ~ field lines reconnect simultaneously 

(|Lazarian and Vishniad [l 9 9 9j) . As a result, the small lo- 
cal reconnection speed V r ' ec loc ps Vas" 1 / 2 is adequate to 
power a GRB if (L/A^s" 1 / 2 > 1, or 



An < L 



2/3 (V\ 1/3 



(45) 



For the Bohm diffusion limit (Eq.[33]), and taking L ~ 
Rg (Eq.pjj]), the condition for triggering a GRB by 
global turbulent reconnection can be written as 



A||<^ /3 r' 1/3 



XT, 



1/3 



B,e 



5° 



1.9 x 10 9 cm 



2/3 



-1/6^1/3 



J w,52 x 2.5 



r 2.5 R U 



1 + a 



1/6 



(46) 



[20fj9l. |2oToIY . Many questions still remain open. 



This requires that the individual reconnection region 
length scale A|| is sufficiently small (e.g. much smaller 

than R/T = 3.3 x 10 12 cmi?i 5 r^, the typical "observ- 
able" scale due to relativistic beaming). This calls for 
multiple collision-induced perturbations to trigger tur- 
bulence as we discuss below. 

3.5. Trigger of a reconnection-turbulence "avalanche" 
and the role of collision 

We envisage internal interactions within the wind (i.e. 
collisions) as the main agent to induce turbulence. Simi- 
lar to the IS picture, numerous collisions can occur within 
the wind. However, most collisions do not lead to a 
strong radiation signature. In the high-cr regime, strong 
shocks can still exist as long as the ram pressure received 
by a magnetize d shell is higher than the rest -frame mag- 
netic pressure (|Zhang and Kobavashil [20051 ) . However, 
if no significant magnetic energy dissipation occurs, the 
total released energy is at most (1 + a)^ 1 of the total 
ejecta energy, even if all the baryonic energy is converted 
to heat and gets radiated away completely. This is the 
main reason that radiation efficiency is low wit hout mag- 
netic dissipation (jZhang and Kobavashill2005l ) . 

In order to give rise to an efficient radiation episode, 
the collision must lead to an abrupt discharge of 
the magnetic field energy, so that the ending a 
value is significantly smaller than the initial one, 
i.e. a en d <IC o-i n i. This must be accompanied by 
an "avalanche" of magnetic reconnection/turbulcncc 
events. A critical condition to trigger such a run- away 
reconnection/turbulence avalanche has to be satisfied. 
A quantitative description of such a critical condition 
is difficult, given the complicated physics involved in 
magnetic reconnection and turbulence. Nonetheless, one 
may envisage the following scenario (FigQ]): 



• Initially, the field lines entrained in the ejecta are 
globally ordered. Since the poloidal field drops 
with radius more rapidly (oc R~ 2 ) than the toroidal 
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field (oc the ordered field lines essentially lie 

in the plane of the ejecta front (i.e. perpendic- 
ular to the direction of motion of the ejecta). In 
the literature, two types of field configurations hav e 
been discussed for GRBs (e.g. ISpruit et al.l (2001): 
the "striped wind" geometry with alternating field 
configuration and the helical "jet-like" structure. 
The former is relevant for pulsar systems in which 
the magnetic axis is misaligned with the rotational 
axis. For GRB systems, unless the central object is 
a rapidly rotating millisecond magnetar, it is more 
likely that the magnetic axis is aligned with the 
rotational axis (e.g. in the black hole - torus sys- 
tems) . This leads to the helical magnetic field con- 
figuration. In the following, we will mostly focus on 
the helical structure, but will mention the striped 
wind geometry when relevant. 

• Internal collisions in an unsteady outflow would al- 
ter the magnetic field configuration. Initially the 
field lines have a large coherent length (Aii ~ Re)- 
Since the most field lines have the same orientation, 
reconnection is very difficult to occur. Even if it oc- 
curs, since the condition to power a GRB (Eqs.[45] 
and [IB]) is not satisfied, the reconnection process 
proceeds very slowly. Without significant magnetic 
dissipation, the collisions are essentially elastic. In 
any case, the ram pressure received by each magne- 
tized shell during the collisions serves as the agent 
to distort the field lines. This is because in real- 
ity, all the fields (density, velocity and magnetic) 
are not uniform, and counter-streaming of the two 
fluids involved in a collision would amplify the per- 
turbations to make the field lines more distorted. 
The scale of coherent length Aii would progressively 
reduce as collisions continue to occur. As a result, 
the parameter s would gradually decrease with ra- 
dius (FigE]). 

• The first reconnection seed requires bringing to- 
gether two sets of field lines with opposite orienta- 
tions. This is easier to achieve for the striped wind 
geometry during the collision processes, but is more 
difficult for the helical geometry. In the majority of 
the locations for the helical configuration, the field 
lines in the two colliding shells have the same orien- 
tation. However, astrophysical systems are not per- 
fectly symmetric systems. One may imagine some 
scenarios to trigger reconnection events at some lo- 
cations. For example, current-driv en kink insta- 
bility may develop in the jet (e.g. iMizuno et al.1 
2009b), which would introduce a slight misalign- 
ment of the magnetic field axes in two consecutive 
"shells" . This would result in a small cross section 
near the magnetic axes that have opposite orien- 
tations in the two shells (Fig|3j) . The first seed 
reconnection would then occur. 

• Once reconnection seed event is triggered, the rapid 
energy ejection from the reconnection layer would 
disturb the nearby ambient plasma and field lines 
to make the region more turbulent. The recon- 
nection layer is also subject to various instabil- 
ities, e.g., the two-stream instability that devel- 
ops as a consequence of reconnection accelera- 



tion. The turbulence stirs the plasma to allow 
more opposite-orientation field lines approach each 
other, so that further reconnection events occur. 
These new reconnection events would eject more 
energy to make the ejecta more turbulent, lead- 
ing to a reconnection-turbulence avalanche. This 
would result in a runaway catastrophic release of 
the stored magnetic field energy. Particles are 
accelerated either directly in the reconnection re- 
gions, or stochastically in the turbulent regions, 
which radiate synchrotron photons that power the 
observed GRB emission. The discharge process 
ceases when a is reduced to a en d < 1 "C cr^. This 
is one ICMART event, which would compose one 
fundamental unit of a GRB (one pulse). Other col- 
lisions that trigger other reconnection-turbulence 
avalanches would give rise to other pulses in the 
GRB. A GRB event is the superposition of many 
individual ICMART events. 

GRB events are relativistic in nature. Although the 
condition (|45|) or (|46|) gives the requirement for pow- 
ering a GRB with global turbulent reconnection. In or- 
der to achieve such a relativistic reconnection/turbulcncc 
avalanche, some physical conditions may be needed. 

The first relevant condition is the strength of the col- 
lisions that distort the field line configurations. We first 
comment on the role of internal shocks in these colli- 
sions. A common misconception is that shocks do not 
exist in the high-cr regime. This is not always true. It 
is relevant for the interaction between a strongly mag- 
netized (high-a) shell and a non-magnetized (a = 0) 
medium (e.g. deceleration of a magnetized shell by the 
ISM). No shock can form in the high-cr shell if the ther- 
mal pressure in the shocked non-magnetized medium 
is smaller t han the magnetic pre s sure i n the magne- 
tized shell (IZhang and Kobavashil l2005t IMizuno et al.1 
l2009aT ) 14 . On the other hand, for the collision between 
two shells with the same a value, even a moderate colli- 
sion would give rise to an excess pressure (ram pressure) 
that exceeds the magnetic pressure. This would lead to a 
pair of shocks passing through both shells. For a given a 
value, the strength of the shock increases with r^ s , and 
can reach more than 50% of the st rength of a a = shock 
once ff., exceeds 3 (Figs 1 & 2 of I Zhang and Kobavashil 
I2005f ). So in the ICMART model essentially all the colli- 
sions are accompanied by internal shocks, but the main 
energy dissipation is not through internal shocks. What 
the shocks do is to disturb the velocity, density and mag- 
netic fields. These distortions are reinforced via turbu- 
lence, which reduce Ay, until an ICMART event is even- 
tually triggered. Without numerical simulations, it is 
hard to quantitatively analyze the role of internal shocks 
in the high-cr regime. Qualitatively, for the same a value, 
a larger relative Lorentz factor between the two shells 

14 The reverse shock forming condition introduced by 
IGiannios et alj {2QC3) is erroneous. These authors claimed that 
there should be no reverse shock when a > 0.3, but the later numer- 
ical simulations performed by the same group suggest s that a wea k 
reverse shock does exist in the a > 0.3 regime (Mirnic a et al.120091). 
The num erical result i s consi stent with Zhang and Kobayashi 
( 2005) and Mizuno et al. ( 2009a) who suggested that the only phys- 
ical condition to define the existence of a reverse shock in a mag- 
netized ejecta is whether the forward shock pressure exceeds the 
magnetic pressure of the ejecta. 
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(Tf s ) tends to introduce a larger perturbation in the or- 
dered magnetic fields. Within the rest frame of shell 1, 
the ram pressure exerted by shell 2 is P ram ,2i — ^\iP2 c2 i 
which would be larger than its own magnetic pressure 
Pbj = B[ 2 /%it = (l/2) / 9 , 1 c 2 cr, when 

r„ >_ <«, 

is satisfied. This is the condition for a strong shock- 
induccd-perturbation to occur. 

A second relevant condition is that the outflow from a 
reconnection event has to be relativistic. This would in- 
duce relativistic turbulence, which is essential to induce 
turbulence-related GRB temporal variability (see more 
in Sect 15. II below). The relativistic extension of Sweet- 
Parker and Petschek scena rios have been discussed by 
iBlackman and Field! (|1994f) and lLvutikov and Uzdenskvl 
(120031) . Based on energy conservation conditions, these 
authors argued that under certain conditions, both the 
inflow and outflow speeds (Vi n and V ou t) can reach 
the relativistic regime. These treatments did not con- 
sider force balance across the reconnection layer, an ef- 
fect becom i ng ve ry important in the relativistic regime. 
iLvubarskvl (|2005|) evaluated this effect and correctly con- 
cluded that Vi n can never achieve the relativistic speed. 
In any case, in long term the condition imposed by en- 
ergy conservation should be satisfied, so that the outflow 
can achieve a rela tivistic speed under certain conditions 
(|Lvubarskvl 120051 M. Lyutikov, 2010, private communi- 
cation). This condition is essential to induce relativistic 
turbulence. The condition for a relativistic outflow may 
be ap proximated as s < a (e.g. iLvutikov and Uzdenskvl 
2003), or, in terms of the Bohm diffusion limit 

A,, < ar' B ^ ~ 9.3xl0 3 cmai7 e ,3i; i 1 5 / 2 2 r2.5^i5 

This condition is more demanding than Eq. (|46l) . which 
means that once relativistic turbulence is excited, the 
turbulent reconnection avalanche can naturally power a 
GRB. 

3.6. Preference of a large emission radius 

The ICMART events occur preferably at large radii 
from the central engine. This is because triggering an 
ICMART event requires that the conditions Eq.(|45|) and 
probably also Eq. (|4"8")) are satisfied, which demands A|| 
being small enough. Early collisions (typically at small 
radii) only serve to distort the field lines and do not di- 
rectly trigger ICMART events. Statistically, shells that 
collide at large radii tend to have more distorted field 
configurations and hence, can more easily trigger IC- 
MART events (FigHJ. 

Observationally, it is not easy to infer the GRB emis- 
sion radius -Rgrb from the MeV data alone. However, by 
combining radiation information from other wavelengths, 
there are three independent, indirect ways to constrain 
-Rgrb- (1) Swift observations suggest that the early X- 
ray afterglo w light curves are domin ated by a steep decay 
component (jTagliaferri et al.ll2005j ). which is best inter- 
preted as the the high-latitude emission o f the prompt 
GRB as the emission is abruptly ceased (jZhang et al.l 
2006) . Within such an interpretation, the duration of the 



tail defines the minimum of Rgrb @j ■ The data suggest 
that .Rgrb is typically large (say, -Rgrb > 10 15 cm) if the 
high-latitude emission int erpretation is valid (|Lyutikovl 
l2006aHKumar et al.ll2007t) . (2) Some GRBs have prompt 
optical detections or upper limits. If the optical emis- 
sion arises from the same region as the MeV emission, 
then the brightness of the optical emission can be used 
to derive a constraint on -Rgrb, based on the argument 
that synchrotron self-absorption should not suppress the 
optical flux. The bright optical prompt emission of 
GRB 080319B requires that the op tical emission region 
has to be greater than 10 1 6 cm (|Racusin et al.l 120081 : 
iKumar and Panaitescull2008l ). so that -Rgrb > 10 16 cm 
is inferred if the MeV emission and optical emission origi- 
nate from the same region (c.f. iFan et al.|[200l IZou et al.l 
2009). A systematic analy sis of other prompt optical 
data (Shcn and Zhang 2009) suggests that data are con- 
sistent with -Rgrb > 10 14 cm. (3) If GeV emission is 
from the same region as the MeV c omponent, as is in - 
ferred from most Fermi/LAT GRBs ([Zhang et al.ll201Q[ ). 
the detection of the high energy photons can be used to 
constr ain both the emission rad ius and the bulk Lorentz 
factor (jGupta and Zhangll2008h . Applying this method 
to GRB 08091 6C, the derived -Rgrb is typically greater 
than 10 15 cm (IZhang and Pe'erll2009T) . We can see that 
the three independent pieces of information generally 
point towards a consistent picture, namely, the GRB 
emission radius -Rgrb is relatively large. This is con- 
sistent with the expectation of the ICMART model. 

4. MERITS AND IMPLICATIONS OF THE ICMART MODEL 

The ICMART model invokes internal collisions as the 
trigger of magnetic reconnection/turbulence avalanche. 
It therefore carries the merits of the internal shock model, 
namely, the emission site is "internal", and the central 
engine activity defines the main variability time scale 
(the pulses in the light curve) 15 . On the other hand, 
it can overcome or alleviate some of the drawbacks of 
the IS model as discussed in §2.1, which we elaborate in 
the following. 

4.1. Radiative efficiency 

Let's assume that two shells [((mi,ri) and (1712,^2)] 
with an initial cr ini collide with each other, and that the 
collision triggers a reconnection/turbulence cascade that 
results in a catastrophic discharge of the magnetic en- 
ergy. This leads to conversion of a significant amount 
of magnetic energy to internal energy of the fluid, and 
hence, to radiation. Let's envisage a picture where the 
two shells merge with a much lower magnetization pa- 
rameter a en d by the end of such an inelastic collision. 
Energy conservation 

(r 2 m 2 +rimi)(l + <Jini) = r m (mi + m 2 +U')(1 + a end ) 

(49) 

and momentum conservation 

(T 2 f32m2+T 1 j3 1 m 1 )(l+o-i n i) = T m /3 m (m 1 +m 2 +U'){l+a en d) 

(50) 

give the same solution for T m [Eq.lJSJ] as the IS model. 
The energy dissipation efficiency, on the other hand, is 

Again this is also the merit of the dissipative photosphere 
model which does not demand internal collisions. 
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much larger, i.e. 



/icmart (r imiC 2+r 2 m 2 c2)(l + a ini ) 
_ 1 T m {m\ + m 2 ) 

1 + CFend (riffli + r 2 m 2 )(l + (Tini) 

^ — (if CTini > 1). (51) 

1 + °~end 

This is ~ 50% if (J en d ~ 1, and can reach 90% if (Tent; ~ 
0.1. This can naturally account for the observed h igh 
radiative efficiency of most GRBs ()Zhang et al.ll2007t ). 

4.2. Particle acceleration and the fast cooling problem 

In an ICMART event, particles can be accelerated 
in the following two ways. First, particles are accel- 
erated directly in the reconnection zone. This hap- 
pens when particles are bounced back and forth be- 
tween the oppositely directed approaching magnetic 
fluxes and gain energy during each gyration. It is sim- 
ilar to shock a cceleration and is an efficient 1st order 
Fermi process (Ide Gouveia dal Pino and Lazarianl [20051 : 
lLazarian and Opherl 12009^ The second acceleration 
mechanism is the stochastic acceleration of particles in 
the turbulent region, which is a 2nd order Fermi process. 

We defer the detailed discussion on particle accelera- 
tion and radiation spectrum of the ICMART model to 
a future work, but qualitatively discuss the difference 
between this model and the shock model. Within the 
standard internal shock model, particles are accelerated 
at the shock fron t via the lst-order Fermi me chanism 
(jSpitkovskvl 120081; iBaring and Summerlinl I2009D . After 
the shock passes through the shell, the downstream hot 
electrons would cool rapidly through synchrotron and in- 
verse Compton emission. Once an electron is acceler- 
ated to a high energy, there is no further heating pro- 
cess, and the electron suffers radiative cooling. Such 
a scenario inevitably introduces the so-called fast cool- 
ing problem as discussed in 32.1.31 unless a much small- 
scale magnetic field is intro duced (jPe'er and Zhang||200fl ; 
lAsano and Ter asawa 20091). 

In the ICMART scenario, electrons do not undergo 
such a "one-shot" acceleration like in the shock scenario. 
Rather, they would be accelerated repetitively. This is 
because the ICMART avalanche happens globally in the 
dissipation region (unlike just near the shock front in the 
shock scenario). There is no distinct separation of the 
particle acceleration region and the cooling region. As a 
result, any particle would undergo multiple accelerations 
in the reconnection regions either directly or stochastic 
scattering of the constantly evolving turbulent field. For 
the lst-order reconnection acceleration scenario, particle 
acceleration time could be much shorter than the cool- 
ing time scale, so that a power-law distributed electron 
population can be accelerated up to an energy at which 
cooling time scale balances the acceleration time scale 
(see Sect 14.31 below). This is similar to the shock sce- 
nario, in which the fast cooling problem remains. On 
the other hand, for the 2nd-order turbulence accelera- 
tion scenario the acceleration time can be long so that 
one may achieve a balance between cooling and heating. 
In the past, this was referred to as the "slow volume heat- 
ing" scenario, which invokes Compton upsc attering off 
emission from a dissipative photosphere (e.g. iThompsonl 
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ous slow heating models is typically sub-relativistic, say, 
lOs-lOOs keV. In the ICMART scenario discussed here, 
the slow heating due to turbulent acceleration would hap- 
pen in the relativistic regime if a significant fraction of 
the magnetic energy is transferred to particle energy via 
the 2nd order stochastic turbulence acceleration. This is 
due to the combination of two effects: a much reduced 
electron number in a high-cr flow (reduced by a factor 
1 + cr) and a much lower magnetic field at a much larger 
emission radius than photosphere (i.e. -R ICMART R P h)- 
To show this, we assume that a fraction £ of comoving 
magnetic energy is consumed through slow heating by 
turbulence acceleration. The average heating rate per 
electron in the emission region is E' heat — £U' B e e / At'n' e , 
where s e is the fraction of the dissipated energy that 
is distributed to electrons. Noticing that the dominant 
cooling mechanism in a high-cr flow is synchrotron ra- 
diation, the average cooling rate can be estimated as 
E' C ooi — (^/3)je a TcU' B . Balancing heating and cooling 
to make E' heat — E' coo[ , one can solve for the critical elec- 
tron energy below which slow heating dominates, i.e. 

xL-^R^(l + a)Y 2 At-r , (52) 

where At' = TAt has been used, and At has been nor- 
malized to 10 ~ 3 s, the minimum variability time scale de- 
tected in GRBs. Equation (f52"j) therefore gives the max- 
imum electron Lorentz factor below which slow heating 
dominates and fast cooling problem is no longer relevant. 
The corresponding synchrotron photon energy is then 

Eh-HT^Ai + z)- 1 ~ 23 keV (&) 
m e c \ Y ) 

x L^T 2 . 5 R 15 m(l + *)i (j^j V2 ACJ (i±£) (53) 

In reality, both acceleration mechanisms would oper- 
ate in the ICMART regions. The lst-order reconnection 
acceleration likely dominates, which would define E p and 
high energy emission of the GRB. If the 2nd-order tur- 
bulence acceleration can transfer a good fraction energy 
to particles so that £ is not -C 1, then Eh is close to E p 
(Eq. [57] discussed in the next section). The fast cooling 
problem may be alleviated. 

4.3. Electron number and E p 

The electron number excess problem of the internal 
shock model ( 32.1.4p is naturally overcome in the IC- 
MART model. This is because the baryon-associated 
electron number is simply (1 + er) -1 of that in the 
baryon-dominated model for the same L w . In princi- 
ple, pair production can increase the lepton number, but 
the large emission site related to the ICMART model 
is way above the pair photosphere dMeszaros et alJl2"002t 
iKobavashi eta l. 2002; Pe'er et al . 2006) so that pair pro- 
duction may not be effective. This makes the Lorentz 
factor of each electron increase by a factor of (1 + a) with 
respect to the baryon-dominated case, naturally raising 
E p to the 7-ray band. 
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To quantitatively address this point, we investigate the 
typical energy of synchrotron radiation in a high-u flow. 
Synchrotron radiation is the dominant mechanism in a 
Poynting-flux-dominated flow, since the photon energy 
density (relevant for inverse Compton scattering) is much 
smaller than the magnetic energy density. The spectral 
peak energy E p in the ICMRT model can also be ex- 
pressed as Eq.([5]). The comoving magnetic field can be 
estimated according to Eq. (j2"0)) . The factor er/(l + a) 
ranges from ~ 1 (for <7j„j 3> 1 in the beginning of an IC- 
MART event) to - 1/2 (for a end ~ 1 at the end of the IC- 
MART event). Assuming that the particle acceleration 
is dominated by the 1st order reconnection acceleration, 
and that the energy release from magnetic dissipation is 
distributed to electrons and protons in the fractions of 
e e and e p with e e + e p = 1 (to be differentiated from e e , 
e p and cb in the shock case), one can write 



L w rj 



AttR 2 cT 2 



^'plprripC 2 + n' e j e m e c 2 = n' e j e m e c 2 /e e . 



( 54 ) 

The comoving electron number density can be calculated 
by Eqs.|HD and CU- Tnis S ives 



7 e = i]e e {l + a)—^ 



Ym e _ 



(55) 



The minimum electron Lorentz factor can be derived by 
7e,m = 0(p)7e- In the ICMART model, the observed 
gamma-ray luminosity can be expressed as 

L 1 = L w rje e . (56) 

one can then write E p as 



E, 



p, ICMART 




Compared with the E p expression in the IS model 
(Eq.g]), the ICMART model can easily make E p in the 
range of several hundreds of keV (as observed) by requir- 
ing a ~ several 10s. 

By reducing the number of electrons, one can also nat- 
urally overcome the self-absorption pr oblem introduced 
by th e electron number excess problem ([Shen and Zhang] 
120091) . 

4.4. Amati/Yonetoku relation 

Comparing Eq.© and Eq. ([57|) . we can see that E p 
of the IS model and ICMART model both have the de- 

1 / 2 

pendence E p oc L 7 . The problem of the IS model (as 
discussed in 8 32.1. 5j) is that Ris has an ad ditional de- 
pendence on T, which is correlated to L 1 (jLiang et al.l 

120101) . This completely destroys the E p cx Ll/ 2 depen- 
dence. In the ICMART model, E p has the dependence of 
oc i 1 / 2 CT 2 i? _1 . It is hard to analytically investigate 

ICMART ° 

1/2 

how a and -R ICMART influence the apparent E p oc L 7 de- 
pendence. The sensitive dependence on a suggests that a 



small variation of a would lead to a significant variation 
of E p . On the other hand, the shallow dependence of a 
on o~o (2/3 power) and R (-1/3 power) suggests that the 
variation of a may not be significant. Also, an outflow 
with a higher a would be more difficult to trigger an IC- 
MART event, since it takes more stringent criteria (e.g. 
a larger V contract and more collisions) to distort the 
field lines (e.g. Eq.gZj). Therefore an ICMART event 
with a higher a tends to happen at a larger R 1CMAWT - 
This would partially compensate the E p scatter due to 
a variation. If the o 2 R~ x factor has a weak depen- 

ICMART " 

dence on L 1 and does not have a large scatter, then one 
may obtain the Amati/Yonetoku correlation. Numerical 
simulations of various processes (e.g. magnetic field dis- 
tortion during collisions, reconnection physics, and the 
trigger condition of ICMART events) are needed to ver- 
ify or disproof such a speculation. 

4.5. Weak photosphere 

Within the ICMART picture, the a value remains 
moderately high at -Rgrb before strong magnetic field 
dissipation happens. This requires that the magnetic 
acceleration and dissipation effects are not prominent 
at the smaller radii. This guarantees a weak photo- 
sphere emission component (energy smaller by a factor of 
(1 + er) -1 with respect to the baryon-dominated model), 
which satisfies the cons traints of the GRB 080916C data 
(jZhang and Pe'erll2009l ) . 

5. PREDICTIONS OF THE ICMART MODEL 

The ICMART model makes a list of predictions that 
may be used to differentiate it from other prompt emis- 
sion models using the observational data. 

5.1. Variability time scales: engine and turbulence 
components 

The most important prediction of the ICMART model 
is that there are two mechanisms that define the ob- 
served GRB variability time scales. Since the emis- 
sion is triggered by internal collisions, the internal wind 
irregularity inevitably leaves an imprint on the light 
curve. This engine-defined variability is similar to that 
of the internal shock model, and is relevant to a broader 
("slow") component of variability. The characteristic 
time scale of this component is defined by the angu- 
lar spreading time at the large emission radius, which 



JT 2 c - 0.3 s R 



ICMART, 15 2.5 



?. Phcnomcno- 



reads R lt 

logically, this corresponds to individual broad "pulses" 
of the light curves, which is visible in some GRBs. On 
the other hand, within each ICMART emission region 
(the broad pulse) , there are many turbulent regions from 
which highly variable emission is released. This corre- 
sponds to a second variability component with very small 
time scales. In the previous GRB models, the variability 
time in all scales are eith er defined by the central en- 
ine (ISari and Piranlll997fl or by relativistic turbulence 
Naravan and Kumarll2009T ). The ICMART model sug- 
gests that the observed GRB variability is the super- 
position of these two components. The two variability 
components may be differentiated via a proper temporal 
analysis of GRB light curves. 

Visually some GRB light curves do show several visi- 
ble slow bumps superposed by rapid variability features 
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(Fig [4]). Po wer density spectrum an alyses using Fourier 
transform (jBeloborodov et all [20 00) . however, did not 
reveal distinct two-component variability distributions. 
Nonetheless, studies using some more sophisticated sta- 
tistical methods have started to reveal the evidence for 
two variability components at least for some GRBs (e.g. 
Shen & Song 2003; Vetere et al. 2006, see also H. Gao 
et al. 2010, in preparation; R. Margutti et al. 2010, in 
preparation) . 

5.2. E p evolution within a pulse 

In the ICMART model, one ICMART event (triggered 
by one collision) corresponds to one pulse in the GRB 
light curve. During each event, since the magnetic en- 
ergy is continuously converted to the particle energy and 
then released as radiation, one may approximately treat 
the plasma as having a continuously decrease with time 
while the lepton number not changed. According to 
Eq. [57] ) , the peak energy in the ICMART model evolves 

as E p oc 'a 2 . Although L 1 increases shallowly during 
the rising phase, the steep dependence on a compensates 
this. Noticing that during the rising phase L 7 only in- 
creases by a factor of a few while a drops 1-2 orders 
of magnitudes, one expects that in general E p decreases 
throughout a pulse (the fundamental radiation unit) in 
GRB prompt emission (Fig]?]), probably with a steepen- 
ing in decay after the light curve peak. Observationally, 
the majority of GR B pulses indeed show a hard-to-soft 
evolution beha vior (jNorris et all 119861: Ifbrd et alJl!995b 
iLu et al1 l2010). which is consistent with the expectation 
of the ICMART model. On the other hand, a small 
fraction of pulses show the so-called "tracking" behavior, 
i.e. the hardness of the spectrum is po sitively correlated 
with the flux variation across the pulse rtGolcnctski i et all 
11983( 1. These pulses cannot be related to an ICMART 
event. It may, however, be interpreted by emission from 
a turbulent eddy, with the rising and decaying compo- 
nents corresponding to entry and exit of the eddy from 
the field of view (r _1 cone). This requires the existence 
of a large scale eddy that can give rise to a broad pulse 
in the GRB light curve. 

5.3. Gamma-ray polarization evolution within a pulse 

Since the ICMART process destroys the ordered global 
magnetic fields in the ejecta, one naturally expects an 
evolution of the linear polarization degree in the j-ray 
emission during each ICMART unit, i.e. a pulse. The 
initial gamma-ray polarization degree may be close to 
the maximum value achievable for synchrotron emission 
in an ordered magne tic field ; i.e. II ~ 50% — 60% 
([Lvutikov et alJl2003t [G rano tl l2003l ). At the end of the 
ICMART event, the ordered field structure would be 
largely destroyed. However, even in the fully turbulent 
limit, the net polarization degree would not become zero, 
since in the strong field regime, turbulence is anisotropic 
with eddies elongated along the local field lines. The av- 
erage polarization in the r _1 -cone field of view would 
not be zero. Without detailed numerical simulations of 
an ICMART event (which requires a dissipative, rela- 
tivistic MHD with the details of reconnection processes 
delineated, and which is beyond the scope of all the cur- 
rent GRB motivated MHD simulations), it is hard to 
quantitatively predict the polarization degree at the end 



of the ICMART event. Nonetheless, a value of n ~ a 
few% — 10% may be reasonable. Since different pulses 
correspond to different ICMART events in our model, we 
expect a significant polarization degree evolution during 
the GRB prompt phase, with II starting with a high value 
at the beginning of each pulse and evolving from high to 
low across the pulse (Fig]?]). For bursts with overlap- 
ping pulses, the II evolution would be more complicated 
since it reflects the superposition of the contributions of 
different pulses. 

Observationally, th ere is no robust gamma-r ay polar- 
ization detection yet. ICoburn and Boggsl (|2003[ ) reported 
n = 80%±20% for GRB 021206 using the RHESSI data. 
However the conclusion cannot be co nfirmed by an in- 
dependent analysis of the s ame data (iRutledge and Foxl 
12004 IWigger et afl 12004( 1. [Willis et all (120051] derived 
lower limits of the polarization degree II > 35% and 
n > 50% for GRB 930131 and GRB 960924, respectively, 
using the BATSE Albedo Polarimetry System (BAPS), 
but could not strongly constrain the degree of polariza- 
tion beyond a systematics-based estimation. Another po- 
larization measurement for the INTEGRAL burst G RB 
041219A (jKalemci et all [20071 iMcGlvnn et all 12007ft is 
also subject to large uncertainty. Nonetheless, these 
studies show the tentative evidence that GRB 7-ray 
emission may be polarized. A true breakthrough could 
be made by one of several proposed 7-ray polarime- 
ters that are suitable to detect polarized 7-ray s from 
transi ent sources, e.g . PO ET ([Hill et alJj2008ft. PO- 
LAR (jLamanna et al.ll2008[ ). and PoGO ([Mizuno et al.l 
12005(1 . For example, for completely ordered magnetic 
fields, POET can measure the polarizations of > 3 0% of 
bursts detected by the mission ([Torna et alJ l2009b). The 
unique time-dependent n evolution pattern predicted for 
the ICMART model may be tested by future observations 
with these detectors. 

5.4. Polarized external reverse shock emission 

Regardless of the composition of the GRB ejecta, 
the outflow is eventually decelerated by a circum- 
burst medium, be it a constant density interstel- 
lar medium (ISM), or a stratified stellar wind 
(iMeszaros and Reeslll997al; ISari et alj|1998t fPai and Lul 
119981: IChevalier and Lil |2000( ). Before entering the 
self-s imilar deceleration regime ([Blandford and McKed 
11976ft . a reverse shock may propagate into the ejecta 
to decelerate it. The composition of the ejecta does 
affect the radiation signatu re from the reverse shock 
(|Zhang and Ko bavashi 2005). In particular, the reverse 
shock emission is expected not to be bright or is com- 
pletely suppressed if a is very small dFan et al.l l2004bl : 
iNakar and Piranl 120041: Din and Fanll2007l) or very large 
(|Zhang and Kobavashi 120051: iMimica et al.l 120091) . A 
bright optical flash may be detected if a < 1, as is dis- 
covered in some G RBs s u ch as GRB 990123 dFan et al.1 
| 2002t IZhang et al.l [20031: iKumar and Panaitescul 120031: 
iGomboc et all 120081 ). At a slightly higher a > 1, the 
reverse shock emission brightness drops significantly. 

In the ICMART model, the final a value after the mag- 
netic dissipation would be close to (maybe slightly below) 
unity, based on an energy equipartition argument. The 
exact a value is hard to predict and is subject to many 
uncertainties. In any case, similar to the discussion of 
prompt 7-ray polarization above, the reverse shock op- 
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tical emission should also be somewhat polarized (say, 
a few percent). If polarimctric observations can be car- 
ried out during the early rapidly decaying phase due to 
the reverse shock emission (say, in the optical band), 
a polarized optical signal should be detected. Obser- 
vations of early optical polarimetry are now regularly 
carried out by some groups, e.g. the Liverpool optical 
polarization observation group. Measurements of early 
optical polarimetry have been made for two GRBs. For 
GRB 060418, which shows a smooth deceleration bump 
that is likely due to the forward s hock emission, an up - 
per limit of II < 8% was set up ([Mundell et alJl2007f >. 
Since the forward shock is the emission of shocked cir- 
cumburst medium, this low level polarization is entirely 
consistent with the theoretical expectation. On the other 
hand, recently the team detecte d a LT = 10% ±2% o ptical 
polarization from GRB 090102 ([Steele et al.ll2009[ ). The 
polarized signal was measured during the rapid decay 
phase of the early optical light curve, which is consistent 
with being the reverse shock emission with a < 1. This 
observational fact is consistent with the expectation of 
the ICMART model. More future early optical polari- 
metric data are needed to further test this prediction of 
the ICMART model. 

5.5. No detectable synchrotron self-Compton spectral 
component during GRB prompt emission 

According to the ICMART model, throughout the 
GRB prompt emission phase, one has a > 1 in the 
emission region. The Compton parameter y is therefore 
<< 1. The magnetic field energy density is much larger 
than the synchrotron photon energy density, so the SSC 
process is greatly suppressed. One therefore expects no 
extra high energy emission features other than extrap- 
olation of the MeV spectrum (Band spectrum) to high 
energy. This is consistent with the observations of the 
majority of Fermi LAT GRBs, including GRB 080916C 
([Zhang et al.ll2010t ). 

6. GRB 080916C 

In this section we discuss GRB 080916C within the 
ICMART model. This is the first bright Fermi/LAT 
burst whose time-dependent broad band spectra are mea- 
sured. The observational results are somewhat surpris- 
ing for theorists. Three theoretically motivated fea- 
tures, namely the pair cutoff feature at high energy, the 
SSC component at high energy, and the quasi-blackbody 
photosphere emission component are all missing. In- 
stead, the time dependent spectra are a set of nearly 
featureless " Band" funct i ons r a nging from ~ 10 k eV to 
- 10 GeV (jAbdo et all [2001 IZhang et all l20Toh . In 
order to interpret this burst within the baryonic mod- 
els, one needs to introduce several spectral components 
that conspire to mimic a Band function covering 6-7 
orders o f magnitude . For example, the model devel- 
oped by iToma et al.l ([2009) invokes SSC in the inter- 
nal shock, upscattered cocoon emission, as well as an- 
other synchrotron emission component to model the ob- 
served Band function spectrum. The neutron- h eatin , 
photosphere model developed by iBeloborodovl ([201 
needs to attribute the very high energy emission to 
a different component (e.g. from the external shock , 
iKumar and Barniol DuranN2009t iGhisellini et al.ll2010D . 
The low energy slope below E p is predicted to be a — 0.4, 



much hard er than the observ ed a ~ — 1 slope in most of 
th e epochs (lAbdo et al.ll2009h . Another model developed 
by Q (|2010a[) applies the standard internal shock syn- 
chrotron emission to account for the emission near E p , 
while invokes a series of residual internal shocks from 
which the IC emission spectra are super posed to account 
for the emission above 100 MeV. Finally, iRazzaqu e et al.l 
(2010) invokes a standard leptonic synchrotron compo- 
nent to interpret the MeV emission, while introduces a 
hadronic proton synchrotron emission component to ac- 
count for emission above 100 MeV. All these models re- 
quire fine tuning of their model parameters to account 
for a very simple broad-band Band spectrum that cov- 
ers 6-7 orders of magnitude in all the time intervals, as 
is r evealed by detailed time-dependent spectral analy- 
ses ([Abdo et all 120091 : IZhang et all l2010f > . Furthermore , 
none of these models have calculated the contribution 
of the photosphere emission (except Beloborodov 2010, 
who discussed the photosphere model itself), but rather 
raised some simple arguments (without calculation) to 
avoid the missing brigh t photosphere problem raise d by 
IZhang and Pe'erT ([2009) . A careful calculation by iFanl 
(2010) showed that the Zhang-Pe'er argument cannot be 
circumvented even under the most favorable condition to 
"hide" the photosphere component. 

The featureless Band function can be straightforwardly 
interpreted within the ICMART model. The lack of the 
three theoretically expected features is understandable 
within the ICMART model. 1. Since the emission radius 
is large for ICMART events ( ^3.6|) . the pair production 
opacity is reduced, so that the rest-frame ~ 70.6 GeV 
photon can escape the emission region. One therefore 
does not expect a pair cutoff feature in the spectrum. 
2. Since the outflow is magnetically dominated in the 
comoving frame (i.e. a > 1), the magnetic energy density 
is naturally much larger than the synchrotron photon 
energy density. The Compton parameter y is <C 1. One 
therefore does not expect a SSC component in the high 
energy regime. 3. The "hot" component only carries a 
luminosity which is (1 + c) -1 times L w . As a result, 
the photosphere emission is much dimmer, which can 
be hidden beneath the no n-thermal emission com ponent 
from the ICMART event ([Zhang and Pe'erll2009D . 

Within the ICMART model, the entire Band function 
spectrum is one emission component, which is powered 
by synchrotron emission of the electrons accelerated in 
the ICMART events. GRB 080916C has an observed E p 
ranging from 400 — 1200 keV in various time intervals 
(|Abdo et all 12009). At z ~ 4.35, its average luminosity 
in the first two main pulses is of the order of ~ (1 — 
2) x 10 53 erg s _1 . From Eq. (|5"T[) . one could estimate the 
required a for this burst 

«7~(38-67)(^ e )- 3 / 4 i^ 
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(58) 

Since (r/e) < 1, the estimate Eq. ([55[) is fully consis- 
tent with the independent constraint a > (15 — 20) 
based on the non-detection of the photosphere compo- 
nent ([Zhang and Pe'erl l2009). One can also estimate the 
maximum synchrotron emission energy in the ICMART 
model. Similar to the shock model, the maximum elec- 
tron Lorentz factor is defined by equating the comov- 
ing acceleration time scale t' acc = K"/ e m e c/ f eB' (where 
k > 1 is a parameter to denote the efficiency of par- 
ticle acceleration) and the electron cooling time scale 
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t' c = Smec/AjeCTTUg (where or is the Thomson scat- 
tering cross section, and U' B = B' 2 /8ir is the comoving 
magnetic field energy density). This gives 

V2 ~ 1.2 x 10 S ^ 2 B'- 1/2 



leM 



6ne 



(59) 



cT'xkB' 

Noticing that the synchrotron spe ctral function F(x) has 
a ma ximum value at x = 0.286 ( Rvbic ki and Lightmanl 
fl979h . the maximum electron synchrotron emission en- 
ergy is 

_3_ ? „ , eB' 

'47T 



: 13 GeV I\«f 



5.35 



(60) 



The expression is simil ar to the shock acceleration case 
(e.g. IWang et all [20091. although k can vary depend- 
ing on the instability growth rate in the ICMART cas- 
cade. In any case, the observed maximum photon en- 
ergy 13.2 GeV can be interpreted given a large enough 
T3 > 1 and a not too large k. A large T is consis- 
tent wit h the constraint deriv ed from the opacity ar- 
gument (jZhang and Pe'erll2009f ). Alternatively, photons 
above 100 MeV may be originated from the external 
shock ( Kumar and Barniol D uran 2009; Ghis ellini et al.l 
l2010t ) 16 . The requirement for k and V is then much less 
demanding. 

Another interesting feature of GRB 080916C is the de- 
layed onset of the LAT band emission. In fact, most 
of th e multi-componen t m odels discussed in the litera- 
ture (|Toma et all 120091: lLl l2010at iRazzaaue et all 120101) 
are motivated to interpret this feature. To us, this fea- 
ture may be straightforwardly interpreted in the follow- 
ing way. The particle acceleration details could be differ- 
ent during the first ICMART event than the later ones. 
Either the electron spectral index is steeper, or there is 
a pair cutoff feature in the LAT band (the latter may be 
related to a smaller T or a smaller emission radius). 

7. CONCLUSIONS AND DISCUSSION 

We have developed a GRB prompt emission model 
in the high-cr regime, namely, the Internal-Collision- 
induccd MAgnetic Reconnection and Turbulence (IC- 
MART) model. This model is motivated by the Fermi ob- 
servations of GRB 080916C (|Abdo et al.ll2009t ). and de- 
veloped upon the idea that GRBs are powered by 3D tur- 
bulence reconnection discussed by lLazarian et al.l (|2003[ ) 
(but in the high-cr regime). This model inherits the 
merits of the internal shock and other models, but may 
overcome several drawbacks of the internal shock model 
(low efficiency, fast cooling, electron number excess, Am- 
ati/Yonetoku relation inconsistency, missing bright pho- 
tosphere component, etc). The basic ingredients of the 
model include the following. 

• The outflow launched from the GRB central engine 
has a high magnetization, i.e. a 3> 1. The a value 

16 This requires a high degree of coincidence since a detailed 
time-dependent spectral analysis with as many time bins as possi- 
ble still reveals a serie s of nearly featurel ess Band-function spectra 
for all the time bins {Zhang et al. 201CJ). In any case, all the ar- 
guments for the ICMART model discussed in this paper are not 
affected even if > 100 MeV photons are of the external shock origin. 



does not decrease significantly before reaching the 
GRB emission radius i?GRB , so that at i?GRB , one 
still has 1 < a < 100. As a result, the photosphere 
emission is not bright, and the observed emission 
is dominated by the non-thermal emission released 
during the ICMART events. 

• The central engine is intermittent, launching an un- 
steady wind with variable luminosity and Lorentz 
factor. The wind interacts internally via collisions. 

• Most collisions at small radii from the central en- 
gine do not result in significant energy dissipation. 
They mainly serve to distort the field lines, making 
the field configuration progressively irregular. 

• At a certain large radius, the condition for turbu- 
lent magnetic field reconnection is satisfied. Re- 
connection events rapidly eject energetic parti- 
cles to the ambient, which further drive turbu- 
lence. This results in a run-away discharge of the 
magnetic field energy in a reconnection/turbulence 
avalanche. This is one ICMART event, which cor- 
responds to one broad GRB pulse. During the mag- 
netic field energy discharge, the a value drops from 
the original value to around unity. 

• A GRB is composed of several ICMART events (i.e. 
broad pulses), each marks a catastrophic event to 
destroy ordered magnetic field lines. The peak en- 
ergy E p is expected to drop from high to low across 
each pulse. The 7-ray polarization degree is also 
expected to drop from ~ 50 — 60% to ~ a few % 
during each pulse. The magnetic field configuration 
at the end of prompt emission is largely random- 
ized, but still has an ordered component. The re- 
verse shock emission is expected to be moderately 
polarized. 

• The GRB light curves should have two variabil- 
ity components, a broad (slow) component related 
to the central engine activity, and a narrow (fast) 
component associated with the relativistic mag- 
netic turbulence. 

This model differs from other magnetic GRB prompt 
emission mo dels proposed in the pas t. The EM model 
proposed by lLvutikov~a nd Blandford] (|2003|) invokes an 
extremely high-cr (a > 10 6 ) at the deceleration radius, 
which might not be realized in nature. The variability 
in this model has no direct connection with the cen- 
tral engine activity, while the evidence of an engine- 
related variability (e.g. those in X-ray flares) is mount- 
ing. A large num ber of GRB magnetic models are in 
the MHD regime (iThompsonl 119941: ISpruit et alj [2001 



Drenkhahn and Spruitl l2002t IVlahakis and Konigll 12003 ; 
Gianniosl 120081 : iKomissarov et al.1 120091 ) . These models 
invoke magnetic dissipation at smaller radii to enhance 
the photosphere emission. At large radii, it is assumed 
that the outflow is no longer Poynting-flux-dominated, so 
that the internal shock model can still operate. On the 
other hand, in the ICMART model it is envisaged that 
rapid reconnection/turbulence cascade only happens un- 
der a certain trigger condition, preferably at a large emis- 
sion radius when the field lines are sufficiently distorted. 
So the main difference between the ICMART model and 
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other MHD models is whether the magnetic energy is re- 
leased abruptly at a large radius or continuously at small 
radii. 

The physics invoked in this model is complicated. In 
this paper, we only limit ourselves to an analytical, qual- 
itative delineation of the general picture of the model. 
Many ingredients of the model, such as magnetic ac- 
celeration in a high-cr flow, collision physics of high-cr 
shells (shocks and magnetic field distortion), reconnec- 
tion physics, particle acceleration (lst-order vs. 2nd- 
order Fermi acceleration) and radiation, are introduced 
based on the best known results from the literature. 
Many speculations are subject to quantitative analyses 
and numerical simulations to verify. Further investiga- 
tions are needed and indeed in plan. 

Within the ICMART model, since the baryon num- 
ber is smaller by a factor of (1 + cr) -1 than the pure 
baryonic model, the expected hadronic radiation is 
also smaller by the same factor. Internal shocks have 
been proposed as th e source of ult rahigh energy cos- 
mic rays (UHECRs, IWaxmanl I1995T) and PeV neutri- 
nos (jWaxman and Bahcalll I1997D . The ICMART pro- 
cesses can in principle also accelerate protons to ultra- 
high energies. However, the UHECR flux from a high-cr 
GRB is smaller by a factor of (1 + cr) -1 than a baryon- 
dominated GRB. If the majority of GRBs have high-cr, 
then GRBs cannot be the dominant contributor to the 
observed UHECRs in the solar neighborhood. Similarly, 
the predicted diffuse PeV neutrino background is also 
lowered by a factor of (1 + cr) _1 . This makes it more 
challenging for km 3 neutrino telescopes (e.g. Icecube) to 
detect these neutrinos from GRBs 17 . 

On the other hand, the GRB composition may be di- 
verse, namely, the a value may vary in a wide range 
among GRBs. For example, another LAT burst GRB 
090902B shows a clear time-evolving blackbody com- 
ponent superposed on a non-thermal power law com- 
ponent (|Rvde et all 120101: IZhang et all I2010D . which is 
almost certainly the baryonic photosphere component. 
This burst is likely orig i nated from a baryonic fireba ll 
pro posed by iPaczvnskil (|1986l ) and iGoodmanl (|1986f) . 
see iPe'er et all (|2010D for a more detailed discussion. 
However, GRB 090902B is a special event. Its special 
spectral f eature is unique wit hin the Fermi LAT GRB 
sample of IZhang et"all (|2010h . Most LAT GRBs have 

17 Our discussion is relevant to traditional high luminosity 
GRBs . The nearby low-luminosity GRBs may be more abun- 
dant l ILiang et al.l [20071; IVirgili et a.1.1 [2009I') . So far there is no 
strong evidence that they are Poynting flux dominated. If they 



clear Band-only time-resolved spectra similar to GRB 
080916C which are good c andidates for the ICMART 
scenario (jZhang et al.ll2010Q . The ICMART model pro- 
posed here therefore can be applied to most GRBs. 

The ICMART scenario developed in this paper may 
be also applied to other astrophysical objects, such 
as active galactic nuclei (AGNs). For example, two 
blazars (Mrk 501 and PKS 2155-30 4) were detected to 
have a 3-5 minute Te V variability (jAlbert et al.l l2007t 
lAharonian et al.ll2007f) . much shorter than the inferred 
light-crossing times at the black hole horizon. The 
detections of the TeV photons also require a Lorentz 
factor much larger than that inferred from the large 
scale jet modeling. The data demand some small- 
scale enhanced emis sion units ([Begelman et al.l 120081 : 
iGiannios et al.l 120091) . Suppose that an intermittent, 
moderately high-cr outflow is launched by the supermas- 
sive black hole central engine from these blazars, IC- 
MART events similar to what are discussed in this paper 
may occur, which would produce small-scale turbulent 
emission units whose size is much smaller than the black 
hole event horizon. These turbulent eddies are relativis- 
tic, and hence, would give an extra Lorentz boost to the 
comoving emission. This would account for the observed 
rapid TeV variability and the apparent large Lorentz fac- 
tor of the TeV emission regions of these two blazars. 
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c speed of light 

c s speed of sound 

d separation between shells in the lab frame, = cSt 

rfmax, d m i n maximum and minimum separation between shells 

e electron charge 

h Planck constant 

tl reduced Planck constant 

k Boltzmann constant or turbulence wave number 

fcj_, fey MHD turbulence wave number in the direction perpendicular /parallel to the magnetic field 

i, /' mean free path of microscopic interactions in the lab frame and comoving frame, respectively 

I' col comoving mean free path of electron Coulomb collision 

m a dimensionless parameter defined as 1 + Ym e /m p 

mi , i«2 mass of shell 1 and shell 2 

m e , rn p electron rest mass, proton rest mass 

n particle number density of the GRB outflow in the lab frame 

n e , n' e electron number density in the lab frame and comoving frame 

n p , n„ proton number density in the lab frame and comoving frame 

n G1 Goldreich- Julian charge number density of the GRB outflow in the lab frame 

p power law spectral index of particles (electrons or protons) 

t 'b e> t 'b p comoving electron gyro-radius and proton gyro-radius 

r co [ strong Coulomb collision radius 

s dimensionless parameter defined in Eq. 1 143 II 

t' acc comoving particle acceleration time scale 

t' c comoving cooling time scale 

t' dyn comoving dynamical time scale 

v' e comoving electron speed 

z redshift 

B, B' magnetic field strength in the lab frame and comoving frame 

B magnetic field vector in the lab frame 

B q critical magnetic field strength 

E observed photon energy 

Eh the critical photon energy below which slow heating effect is important 

Em maximum observed photon energy 

E electric field vector in the lab frame 

E p observed spectral peak energy of GRB 

E lt i ao isotropic gamma-ray energy of GRB 

E(k) turbulence energy per unit wave number at the wave number k 

E(k±), E(k\\) turbulence energy per unit wave number at the wave number k±, and fcy 

E(8) GRB effective "isotropic" energy at an angle 9 from the jet axis 

B'heat' ^'cool comoving average heating and cooling rate of electrons 

F(x) a function to denote synchrotron emission spectrum of a single particle 

Ft, matter flux: sum of baryonic and leptonic fluxes. Usually baryonic flux dominated. 

Fp Poynting flux 

J current density vector 

L characteristic length scale of the flow 

C length of reconnection layer 

L w the total (including kinetic and magnetic) isotropic luminosity of the GRB ejecta (wind) 

L 1 the isotropic 7-ray luminosity of the GRB 

N(E) number of photons in the energy bin (E, E + dE) 

7V(7 e ) number of electrons in the Lorentz factor bin (-y e , 7 e + d-y e ) 

P pressure 

Pgas gas pressure 

Pmag magnetic field pressure 

Pram, 21 ram pressure exerted to shell 1 by shell 2 

R radius from the central engine 

Rdcc GRB ejecta deceleration radius 

R e Reynold's number 

i?GRB radius of GRB prompt emission from the central engine 

iJjCMART radius of the ICMART events 

RlS internal shock radius defined in Eq. (2j 

R m magnetic Reynold's number 

-R MHD radius from central engine where the MHD condition is broken 

Rg jet cross section radius at the radius R 

S Lundquist number 

T characteristic temperature of plasma 

T e , T p characteristic electron temperature, proton temperature 

T e , c critical electron temperature for collisional/collisionless regime separation 

U 1 total internal energy in the internal shock or ICMART event 

U'g comoving magnetic field energy density 

Up h comoving photon energy density 

V velocity vector of the GRB fluid 

Vj n incoming speed of the magnetic field lines in Sweet-Parker reconnection 

Vout the eventual outgoing speed in a reconnection event 

V\ Alfven speed in general 

VI comoving Alfven speed 

V£ NR comoving Alfven speed in the non-relativistic regime 
K'cc.loc' K'cc, global comoving local and global reconnection speed 

x argument of function F(x) 
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Y number of leptons associated with each proton 

y The Compton parameter 

a GRB photon spectral index below E p 

P dimensionless speed V/ c or ratio between gas pressure and magnetic field pressure 

f dimensionless speed of the fast shell 

m dimensionless speed of the merged shell 

Pa dimensionless speed of the slow shell 

7^1 Alfven Lorentz factor in general 

7^ comoving Alfven Lorentz factor 

7e comoving electron Lorentz factor 

7e, c comoving electron Lorentz factor at cooling break 

7e comoving mean electron Lorentz factor 

7m i 7m minimum and maximum Lorentz factor of a power-law distributed proton or electron population 

7e,m minimum comoving electron Lorentz factor in an injected power law energy spectrum 

7 e m maximum comoving electron Lorentz factor in an injected power law energy spectrum 

7e lP comoving electron Lorentz factor that contributes to E p in the synchrotron radiation model 

7 P comoving proton Lorentz factor 

7 P comoving mean proton Lorentz factor 

7,; n Lorentz factor of incoming magnetic field lines for relativistic reconnection 

r GRB ejecta bulk Lorentz factor 

To the ejecta bulk Lorentz factor at the central engine 

Ti , T?2 Lorentz factor of shell 1 and shell 2 

F21 relative Lorentz factor between shell 2 and shell 1 

T43 relative Lorentz factor between regions 4 (unshocked trailing shell) and 3 (shocked trailing shell) 

r i Lorentz factor of the fast shell 

r m Lorentz factor of the merged shell 

r s Lorentz factor of the slow shell 

Ty s relative Lorentz factor between the fast and slow shells 

Y ul i relative Lorentz factor between upstream and downstream 

r max maximum Lorentz factor in the GRB ejecta 

r m i n minimum Lorentz factor in the GRB ejecta 

Ttot achievable Lorentz factor of high-cr shell for total conversion of Poynting energy to kinetic energy 

8 thickness of reconnection layer 

&' e , &' p comoving electron and proton plasma skin depth 

8t duration between the end of ejecting a leading shell and the beginning of ejecting a trailing shell 

<5<max maximum St in the ejecta 

<5t m ; n minimum 8t in the ejecta 

8V relative velocity of the flow 

A shell width in the lab frame 

A' shell width in the comoving frame 

A f width of the fast shell 

A s width of the slow shell 

A ma x maximum shell width in the GRB ejecta, = cAt ma x 

A m i n minimum shell width in the GRB ejecta, = cAt m i n 

At duration of central engine activity for each mini-shell in the ejecta 

At' comoving time scale for global magnetic dissipation within a shell with comoving width A' 

Af. max maximum At in the ejecta 

At m i n minimum At in the ejecta 

eg fraction of internal energy that is distributed to magnetic fields in internal shocks 

e e , e p fraction of internal energy that is distributed to electrons and protons in internal shocks 

e e , e p fraction of dissipated magnetic energy distributed to electrons and protons in an ICMART event 

r\ magnetic diffusion coefficient or energy dissipation efficiency in general 

r] 13 energy dissipation efficiency of an internal shock 

^icmart energy dissipation efficiency of an ICMART event 

6 angle from the GRB jet axis 

6j GRB jet opening angle 

k a parameter to denote efficiency of particle acceleration 

\g coherence length of a random magnetic field 

A || local reconnection length 

v kinematic viscosity 

£ fraction of comoving magnetic field energy that is dissipated through slow heating mechanism 

El linear polarization degree 

p mass density in the lab frame 

p' mass density in the comoving frame 

p' ± , pj comoving mass density of the shell 1 and shell 2 

cr magnetization parameter as defined in Eq. JTJ 

(jo initial cr value at the GRB central engine 

cr c critical cr value to separate the sub-Alfven regime from the super-Alfven regime 

<r en[ j final cr value after an ICMART event 

<Tj n j initial a value before an ICMART event 

cry Thomson scattering cross section 
comoving Coulomb collision time 
non-relativistic comoving Coulomb collision time 
relativistic comoving Coulomb collision time 
viscous diffusion time 

r^i f magnetic resistive diffusion time 

Tr flow time scale 

<f>(p) a function of p to connect minimum particle energy with the mean energy 
u)' B e , u)' B comoving electron and proton gyro-frequency 
u)' p e , u' PtP comoving electron and proton plasma frequency 



col 
T col,NR 

t' , p 

cot.H 
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Notes added in proof: After acceptance of our paper, we were notified by Jon McKinney about an alternative 
mechanism to trigger fast reconnection at large radii of GRBs (J. C. McKinney & D. A. Uzdensky, MNRAS, submitted, 
larXiv: 1011 .l903j) . 
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Distance Scales in the ICMART Model 




central engine photosphere early collisions ICMART region External shock 

R~10 7 cm R~ 10 11 - 10 12 cm R ~ 10 13 - 10 14 cm R - 1 15 - 1 16 cm R~10 17 cm 

o = a »1 oso a~ 1-100 a ini ~ 1-100 osl 

a end < 1 



Fig. 1. — A cartoon picture of the ICMART model. The typical distances and a values of various events are marked. 
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(a) Initial collisions only distort magnetic fields 




(b) Finally a collision results in an ICMART event 



Fig. 2. — Examples of internal collisions that mainly distort magnetic fields and result in catastrophic discharge of magnetic energy in 
an ICMART event. 
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Fig. 4. — An example of GRB light curve that shows two variability time scales. The light curve of GRB 050117 is extracted from 
the UNLV GRB group website http://grb.physics.unlv.edu/~xrt/xrtweb/050117/050117.html The predictions of decreasing gamma-ray 
polarization degree n and the spectral peak energy hJ p within individual pulses are indicatively presented. Detailed decaying functions 
would be different depending on the details of evolution of magnetic field configuration, a value, as well as balance between heating and 
cooling of electrons. The general decreasing trend is robust. 



